Astronomy & Astrophysics manuscript no. hisraer400dir proved 


©ESO 2012 


October 9, 2012 





The 400d Galaxy Cluster Survey weak lensing programme: 
II: Weak lensing study of seven clusters with MMT/Megacam * 

Holger Israel 1 , Thomas Erben 1 , Thomas H. Reiprich 1 , Alexey Vikhlinin 2 , Craig L. Sarazin 3 , and Peter Schneider 1 

1 Argelander-Institut fiir Astronomie, Auf dem Hiigel 71, 53121 Bonn, Germany 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

3 Department of Astronomy, University of Virginia, 530 McCormick Road, Charlottesville, VA 22904, USA 

Accepted for publication in Astronomy & Astrophysics 

ABSTRACT 

Context. Evolution in the mass function of galaxy clusters sensitively traces both the expansion history of the Universe and cos- 
mological structure formation. Robust cluster mass determinations are a key ingredient for a reliable measurement of this evolution, 
especially at high redshift. Weak gravitational lensing is a promising tool for, on average, unbiased mass estimates. 
Aims. This weak lensing project aims at measuring reliable weak lensing masses for a complete X-ray selected sample of 36 high 
redshift (0.35 < z < 0.9) clusters. The goal of this paper is to demonstrate the robustness of the methodology against commonly 
encountered problems, including pure instrumental effects, the presence of bright (8-9 mag) stars close to the cluster centre, ground 
based measurements of high-z (z~0.8) clusters, and the presence of massive unrelated structures along the line-sight. 
Methods. We select a subsample of seven clusters observed with MMT/Megacam. Instrumental effects are checked in detail by cross- 
comparison with an archival CFHT/MegaCam observation. We derive mass estimates for seven clusters by modelling the tangential 
shear with an NFW profile, in two cases with multiple components to account for projected structures in the line-of-sight. 
Results. We firmly detect lensing signals from all seven clusters at more than 3.5cr and determine their masses, ranging from 10 14 M 
to 10 15 M , despite the presence of nearby bright stars. We retrieve the lensing signal of more than one cluster in the CL 1701+6414 
field, while apparently observing CL 1701+6414 through a massive foreground filament. We also find a multi-peaked shear signal in 
CL 1641+4001. Shear structures measured in the MMT and CFHT images of CL 1701+6414 are highly correlated. 
Conclusions. We confirm the capability of MMT/Megacam to infer weak lensing masses from high-z clusters, demonstrated by the 
high level of consistency between MMT and CFHT results for CL 1701+6414. This shows that, when a sophisticated analysis is 
applied, instrumental effects are well under control. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individuals: CL0159+0030, CL0230+1836, CL 0809+2811, 
CL 1357+6232, CL 1416+4446, CL 1641+4001, CL 1701+6414 - Cosmology: observations - Gravitational lensing - X-rays: galax- 
ies: clusters 



1. Introduction 

Cosmological observables probing different physics are found 
to agree within their uncertainties on a ACDM cosmologi- 
cal model dominated by Dark Energy and Dark Matter (e.g., 
|Kowalski et al.|2008||Schrabback et al.|20lQ||Larson et al.|201 1) . 
Investigating the unknown physical nature of Dark Energy ranks 
among the foremost questions for cosmologists. In particular, the 
presence or absence of evolution in Dark Energy density is ex- 
pressed by the equation-of-state parameter wde- State-of-the-art 
measurements (e.g., Vikhlinin et al. 2009bj Mantz et aE]|2010 



Kom atsu et al.||2011| l are consistent with wde = — 1, and hence 
with Dark Energy being an Einsteinian cosmological constant 
(e.g. iBlanchard 2010). The tightest constraints on w DE can be 
achieved only by combining results from different probes show- 
ing complementary dependencies on cosmological parameters. 

The evolution of galaxy clusters is understood to be de- 
termined by cosmological parameters through both cosmic ex- 

* Observations reported here were obtained at the MMT Observatory, 
a joint facility of the Smithsonian Institution and the University of 
Arizona. 

** Reduced and coadded MMT image files are available in elec- 
tronic form at the CDS via anonymous ftp to cdsarc . u-strasbg . fr 
(130.79.128.5). 



pansion and hierarchical structure formation. Thus, clusters pro- 
vide information to constrain cosmological parameters that is 
complementary to other tests (using e.g., the cosmic microwave 
background, type la supernovae, or baryonic acoustic oscilla- 
tions). Information about the expansion history of the Universe 
and structure formation is encoded in the cluster mass function 
n(M,z); different cosmological models imply a different cluster 
mass function at high z, compared to the local mass function 
(e.g. |Eke et al.||T996| |Rosati et al.||2002"t |Voit||2005) [Schuecker 



2005 



Press & Schechter ( 



Reiprich||2006|l. Re fining the original analytic model of 
974)), various fitting formulae have been 



developed based on numerical simulations (see Pillepich et al 
|2010| for an overview). The measured cluster mass function not 
only gives strong evidence for the existence of Dark Energy (e.g. 
Vikh linin et al. 2009b ) but also adds valuable information to the 
joint constraints of cosmological parameters (for a recent review, 
see |Allen etal.|2011| ). 

Galaxy cluster cosmology relies on the accurate determina- 
tion of cluster masses and a thorough understanding of the dif- 
ferent mass proxies in use. Cluster masses are most commonly 
inferred from a variety of X-ray observables (X-ray luminosity 
L x , gas mass M gas , the quantity Y x - 7xM gas ) or gravitational 
lensing (both weak and strong), but as well using the Sunyaev- 
Zel'dovich effect or the motions of member galaxies. For distant 
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galaxy clusters which have low masses owing to the early epoch 
of structure formation they represent, the small number of avail- 
able photons prohibits a detailed (spectral) X-ray analyses avail- 
able for local clusters. Nevertheless, weighing a large number of 
high-redshift clusters will yield the best constraints on cosmo- 
logical parameters. This strategy will be adopted by the upcom- 
ing generation of cluster surveys, e.g. using eROSITA ( |Predehl| 
|et al . |20 10} |Pillepich et al.|201 1 1. Hence, scaling relations con- 
necting quantities like Lx or M gas with the total cluster mass will 
continue to play important roles and need to be understood and 
calibrated thoroughly. 

Rooted in the thermodynamics of the intracluster medium 
(ICM), X-ray methods rely on assumptions of hydrostatic equi 
librium, elemental composition, and, to a large extent, spheric 



ity of the cluster's gravitational potential well (Sarazin 1988 
Bohringer & Werner 2010). Weak gravitational lensing (WL; 



e.g. Schneider 2006] offers an alternative avenue for determin- 
ing cluster masses, which is completely independent of these as- 
sumptions, directly mapping the projected mass distribution of 
matter, Dark and luminous. 

Merging clusters deviate strongly from thermal and hydro- 
static equilibrium, with a significant amount of the internal en- 
ergy being present as kinetic energy of bulk motions or turbulent 
processes, e.g. merger shocks. Thus, the expectation from nu- 
merical simulations is that X-ray masses for merging systems 
might be biased after a significant merger, with a relaxation 
timescale of 0(1 Gyr). Simulations agree with the expectation 
from the hierarchical structure formation picture that mergers 
are more frequent at higher redshifts than in the local Universe 



(e.g. Cohn & White|2005). There is no consensus between differ- 



ent simulations yet which of several suggested physical effects 
dominates after the phase in which a disturbed morphology can 
be seen (e. g., |Kravtsov etaL|2006| |Nagai et al.]|2007a| |Stanek 



|et al.|20i0| ). Most importantly, bulk motions induce non-thermal 
pressure, providing support for the gas against gravity, thus pos- 
sibly leading the hydrostatic mass to underestimate the true mass 



by 5-20% even in relaxed clu sters (Rasia et al. 2006 Nagai et al. 
2007b[|Meneghetti et al.|20T0 



Therefore, studying scaling relations of X-ray observables 
with weak lensing masses has become an important ingredient 
in refining cluster masses from X-ray observations (e.g., Zhang 



et al.|2008||2010||Meneghetti et al.|2010) . Relative uncertainties 



of the individual WL cluster masses are higher than those from 
X-rays, largely due to intrinsic shape noise. But the power of 
weak lensing comes through the statistical analysis of M w \/Mx 
for the whole sample, under the assumption that WL mass es- 
timates are, on average, unbiased. This means, while WL mass 
estimates for individual clusters are subject to an error due to 
the projection of filaments or voids along the line-of-sight, the 
stochastic nature of these errors makes them cancel out when 
averaging over a well-defined cluster sample. Statistical com- 
parisons to X-ray masses (e.g., Meneghetti et al. 2010 ) help us to 
investigate WL systematic uncertainties, i.e. tri axiality (|Corless| 
|& King|2009| and projection of unrelated LSS ( |Hoekstra|2003| ), 
to which X-ray observables are far less sensitive. 

This article presents the second part of a series on weak 
lensing analyses following up the 400 Square Degree Galaxy 
Cluster (400d) Survey, initiated in |Israel et"aT1 ( |2010| hereafter 
Paper I). The 400d Survey presents a flux-limited sample of 
galaxy clusters detected serendipitously in a re-analysis of all 
suitable ROSAT PS PC pointings ([Burenin et al.|2007| ). From the 
resulting catalogue, Vikhlini n et al.| ( |2009a| V09) drew the cos- 
mological subsample of 36 X-ray luminous and distant clusters, 
for which high-quality Chandra X-ray observations were ob- 



tained and analysed. The CHANDRA-based cluster mass function 
resulting from the Chandra Cluster Cosmology Project was pub- 
lished by V09, for the complete redshift range of 0.35 < z < 0.90 
spanned by the clusters in the cosmological subsample, as well 
as divided into three redshift bins. Building on this mass func- 
tion, Vikhlinin et al. ( 2009b ) constrained cosmological parame- 



ters, in particular wde- 

Determining accurate weak lensing masses for the distant 
clusters in the 400d cosmological subsample opens the way to 
observationally test the assumptions Vikhlinin et al. ( |2009a|b[ ) 
make for the scaling relations and their evolution. Put briefly, 
the WL follow-up of the 400d cosmological sample clusters pro- 
vides us with a control experiment for the resulting X-ray mass 
function. With 36 clusters, the 400d WL sample ranks among 
the largest complete high-z WL samples. 

In Paper I, we presented the results of our feasibility study, 
performing a detailed lensing and multi-method analysis of 
CL 0030+261 8. In particular, we showed the Megacam instru- 
ment to be well suited for WL studies. As the next step of the 
project, we investigate seven further clusters from our sample, 
all of which were also observed with Megacam at MMT The re- 
sulting WL mass determination and the status after 8 out of 36 
clusters have been analysed are the subjects of this paper. 

We consistently assume a ACDM cosmology specified by 
the dimensionless Hubble parameter h = 0.72 and matter and 
Dark Energy density parameters of Q m = 0.30 and £2a = 0.70. 

This paper is organised as follows: After giving a short 
overview on our data set and its reduction in Sect. [2] we give 
salient details of the WL analysis methods we used in Sect. [3] 
In Sect. |4] we take a closer look at two clusters which show 
a more complicated shear morphology and devise a two-cluster 
shear model. Comparing our MMT results to a CFHT weak lens- 
ing analysis of one of our clusters, we once more prove MMT 
weak lensing to be reliable and provide an external calibration 
(Sect. [5]l. In Sect. [6] we provide details of the error analysis for 
our main results, the cluster masses, which are then discussed in 
Sect. [7] Section|8]presents our summary and conclusion. 



2. Methodology 

2.1. MMT/Meqacam data for the 400d WL survey 

Table [T] summarises the observations of the eight 8 > 0° galaxy 
clusters with right ascension^] h < a < 8 h 30 m and 13 h 30 m < 
a < 24 h for which MMT/Megacam observations in the lensing 
(r'-) band have been completed. As CL 0030+26 18 was studied 
in detail in Paper I, this work focusses on the remaining seven 
clusters. Following the observation strategy described in Paper I, 
with nominal exposures of r nom = (7500 s, 6000 s, 4500 s) in 
(g'r'i'), these seven clusters were observed in the four out of 
five MMT observing runs performed for the 400d WL survey in 
which weather conditions permitted usable observations during 
at least parts of the scheduled time. 

Megacam ( |McLeod et al.||2000| ), then located at the Fred 
Lawrence Whipple Observatory's 6.5 m MMT telescope, is a 
high-resolution (0.08" px~'), wide-field (~ 24' x 24' field-of- 
view) camera, consisting of a 4 x 9 CCD mosaic. 

The four "winter" clusters CL 0030+2618, CL 0159+0030, 
CL 0230+ 1836, and CL 0809+28 11 have completed observa- 
tions in the g'r'i' filters, while due to scheduling constraints, 
only the r'-imaging could be completed for the "summer" 
clusters CL 1357+6232, CL 1416+4446, CL 1641+4001, and 



Their J2000 coordinates are given by the designations in Table[T| 
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Table 1. Specifications of the data sets for all eight clusters analysed so far. For each cluster and filter, the dates of observation, total 
Megacam exposure time r™ p , usable final exposure time T^ p , seeing, and 5cr limiting magnitude (Eq. 2 in Paper I) for the final 
image stack are given. The last column refers to a direct (D) or indirect (I) photometric calibration (PhC, Sect. 2.2 1. CL 0030+2618 
is the cluster analysed in Paper I. 



Cluster 


Filter 


Observation Dates 


r™ \s\ 


CM 


Seeing 


"'lira 


PhC 


CL 0030+26 18 


r' 


2004-10-06/07 


15300 


6600 


07 82 


25.9 


I 




g' 


2005-10-30/31,11-01 


9150 


7950 


a/87 


26.8 


D 




i' 


2005-10-31 


6000 


5700 


1'.'03 


25.1 


D 


CL0159+0030 


r> 


2005-10-30/31,11-01 


9900 


3600 


a/85 


25.7 


D 




g' 


2005-11-01, 


6000 


4800 


1705 


27.7 


D 




i' 


2005-10-31,11-01 


8100 


5700 


1714 


25.0 


D 


CL 0230+ 1836 


r> 


2004-10-06/07; 2005-11-08 


9600 


2700 


a/68 


25.1 


I 




g' 


2005-11-08 


6000 


4200 


a/ 80 


27.2 


I 




V 


2005-10-31,11-01/08 


9600 


3600 


a/98 


24.7 


D 


CL 0809+28 11 


r' 


2005-11-08; 2008-01-09 


9300 


3000 


a/72 


25.4 


D 




g 1 


2005-10-31/11-08 


6000 


3600 


1704 


26.3 


D 




i' 


2005-10-31/11-01 


7500 


5700 


a/82 


26.1 


D 


CL 1357+6232 


r' 


2005-06-07 


7200 


2700 


a/90 


25.4 


D 


CL 1416+4446 


r' 


2005-06-08 


7500 


4200 


a/ 8i 


25.8 


D 


CL 1641+4001 


r' 


2005-06-07 


8100 


6900 


a/9i 


26.0 


D 


CL 1701+6414 


r' 


2005-06-08 


7500 


6000 


a/89 


25.8 


D 



CL 1701+6414. Therefore, a different strategy has to be adopted 
for parts of the data reduction (Sect. 2.3 i and the background 
source selection (Sect. |2.4| i for these single-band clusters com- 
pared to three-band clusters. 

As indicated in Table [T] some clusters were observed in 
the same filter in more than one observing run. Using the data 
reduction described in Sect. |2.2| we produced coadded (stacked) 
images, for which net exposure times 7*1, seeing, 5<x limiting 



very bright stars is not removed by THELI, but regions in the 
image in which source counts deviate significantly from the 
mean are masked as unreliable using the algorithm described 



magnitudes, and photometric calibration method (Sect. 2.2 1 are 
given in the four last columns of Table [T] 

The most striking fact to note are the drastic reductions 
from the initial raw data exposure times T™ p to the used 
in the coadded images. In a number of cases, the required see- 
ing in the coadded image of < 1" in the lensing band and < 172 
in the other bands could only be achieved by removing images 
such that 7™, < T nom . As this inevitably reduces the limiting 
magnitude (Eq. 2 in Paper I), the final stacks represent a com- 
promise between seeing and depth, aiming at an optimal WL 
signal. Similarly, compromises had to be made between main- 
taining a low level of anisotropy (Sect. |2.3| l in the point spread 
function (PSF) and limiting magnitude. The ramifications of the 
heterogeneous data quality and the - in some cases - shallow 
exposure times, for which the good overall r'-band seeing could 
be obtained, will be addressed at several occasions in this article. 



2.2. Data reduction and calibration 

The data reduction for the 400d WL survey has been described 
in detail in Paper I. Therefore, we give only a brief recapitulation 
here and refer the interested reader to Paper I. 

The first stage, including all tasks of elementary data reduc- 
tion (de-biasing, flatfielding, de-fringing, construction of weight 
images, astrometry, relative photometry, and coaddition) are per- 
formed using the THELI pipeline for optical data reduction intro- 
duced by Erben et al. ( 2005) 1 and adapted to MMT/Megacam in 
Paper I. Generally, we achieve a high level of homogeneity in 
the noise backgrounds of coadded images, with our pipeline 
effectively correcting the position-dependent transmissivity of 
Megacam filters (Appendix [Al. The additive stray-light from 



by Dietrich et al. (2007). The mask images we produce in the 



second stage of data reduction for each coadded image also con- 
tain masks for saturated stars (cf. Paper I) and a few manually 
inserted masks for, e.g., asteroid trails^] The first seven panels 
of Fig. [TJpresent the central regions of our clusters as observed 
with MMT/megacam. For the three-band clusters, we prepared 
pseudo-colour images using the g'r'i' coadded images. 

Applying the method of Hildebra ndt et al.| ( |2006[ ), we com- 
puted absolute photometric zeropoints for our coadded images. 
As photometric reference for this calibration, the Sloan Digital 
Sky Survey Data Release Six (SDSS DR6, |Adelman-McCarthy 



|et al.|2008| ) was employed, with which six of our clusters over- 
lap. This direct calibration also yields zeropoints for fields out- 
side the SDSS footprint observed in the same filter in the same 
photometric night, as for the ;'-band of CL 0230+1836 (Table[T}. 
The remaining observations were done in nights in which no 
cluster with SDSS overlap was observed under photometric con- 
ditions. To these data, labelled with "I" in Table[T] we applied an 
indirect calibration described in Paper I, basically a rudimentary 
but effective stellar locus regression ( |High et aL||2009| l. Details 
concerning the results and accuracy of the photometric calibra- 
tion can be found in Appendix |A| 



2.3. From images to shape catalogues 

A detailed description of how we distill from a coadded image 
a galaxy shape catalogue, containing positions, ellipticity mea- 
surements, and photometric data for sources that can be con- 
sidered galaxies can be found in Paper I. For the single-band 
clusters, we use straight-forward calls to SExtractor ( |Bertin| 
|& Arnouts||1996| l. For the three-band clusters, the convolu- 
tion of images to the seeing in the poorest band and calls to 
SExtractor in double-detection mode, using the (unconvolved) 



2 The more common satellite streaks are masked already during the 
basic data reduction, prior to coaddition. 
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Fig. 1. Clusters discussed in this paper. We show pseudo-colour images for the cases where colour is available, using the MMT 
g'r'i' bands. For the CL 1701+6414 field, we also show pseudo-colour images using the CFHT g'r'i' bands, both for CL 1701+6414 
and A 2246. We choose the rosat cluster coordinates as centre of the images. Note the variable background due to bright stars near 
CL 0159+0030 and CL 0809+28 1 1 . 



r'-band image as the "detection image", as described in Paper I 
are being performed. 

We apply the "TS" shear measurement pipeline (Heymans 
et al.|2006[|Schrabback et al.|2007[ |Hartlap et al.|2009), an"m> 
plementation of the KSB+ algorithm ( |Kaiser et al.|1995| [Erben 
et al.|200l) , which determines moments of the brightness distri- 
bution for each source and corrects for the convolution with an 
anisotropic PSF The PSF anisotropy is traced by measuring the 
brightness distribution of sources identified as stars in a plot of 
their magnitude ^ UT0 against the half-light radius §. The values 
we used to define the boundaries of the stellar locus are given in 



Table [2] Only sources in the KSB catalogue, consisting of detec- 
tions with a viable measurement of §, are further considered. 

Consistent with our Paper I findings, MMT/Megacam ex- 
hibits a smooth, albeit variable pattern of PSF anisotropy which 
can be modelled by a low-order (2 < J an ; < 5, see Table |2|i 
polynomial in image coordinates such that the residual PSF 
anisotropy has a practically vanishing mean value and a disper- 
sion 0.005 < <r(e anl ) < 0.010 in the r'-band image stacks. In terms 
of the uncorrected PSF anisotropy, however, there are consider- 
able differences in the input images for different cluster fields. 
Excessive PSF anisotropy observed in several input frames - 
which thus had to be removed from the coadded images - can 



4 



Holger Israel et al.: The 400d weak lensing survey II 



and r[* ax delineate the stellar locus. 



Table 2. Definitions of the galaxy shape and lensing catalogues. The parameters &* min , #* nax , ' 

The galaxy shape catalogue considers sources i?> ??J nax for r^ jn < r' Am0 < r max an ^ ^ > ^rrJn f° r r AUTO- r max (Cols. (1) to (5)). The 
number densities «ksb and n ga i of sources in the KSB and galaxy shape catalogues are shown in columns (6) and (7); followed 
by the degree c/ an ; of the polynomial for PSF anisotropy correction. Columns (9) to (13) present the photometric cuts mbright and 
fifaint, defining the lensing catalogue with a source density n\ c . Using a filter scale of 6 oM in Eq. ([TJ, we find a maximum shear peak 
significance S max . CL 0030+26 18 (Paper I) is included for completeness. 



Cluster field 


mm 


r 1 

mm 


ft* 

mux 


r". 

mm 


r'* 

max 


"KSB 






^bright 






„opl 
tf out 


c 

u max 




[px] 


[px] 


[px] 


[mag] 


[mag] 


[arcmirT 2 ] 


[arcmirT"] 




[mag] 


[mag] 


[arcmirT - ] 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


CL 0030+26 18 


2.55 


2.80 


2.95 


16.75 


22.5 


44.79 


22.64 


5 


20.0 


22.5 


21.28 


15:5 


5.84 


CL 0159+0030 


2.65 


2.95 


3.10 


17.5 


22.5 


37.27 


18.91 


3 


20.0 


24.6 


11.58 


10:5 


4.12 


CL 0230+1836 


2.05 


2.33 


2.45 


16.75 


22.5 


29.28 


13.22 


4 


21.5 


23.0 


11.04 


8:5 


3.64 


CL 0809+28 11 


2.35 


2.66 


2.80 


17.0 


22.5 


33.39 


16.56 


2 


22.5 


24.8 


10.31 


13:5 


5.39 


CL 1357+6232 


2.85 


3.18 


3.35 


17.5 


22.25 


29.52 


14.24 


5 




18.5 


14.23 


10:17 


4.47 


CL 1416+4446 


2.55 


2.80 


2.95 


16.75 


22.5 


42.50 


21.51 


5 




18.5 


21.50 


4183 


4.25 


CL 1641+4001 


2.90 


3.135 


3.30 


16.75 


22.5 


39.95 


19.58 


3 




22.7 


17.65 


i6:o f 


4.12 


CL 1701+6414 


2.60 


2.95 


3.10 


16.75 


22.5 


40.71 


20.72 


4 




21.9 


19.47 


15:5' 


3.75 



t No significant decrease of S(9 ollt ) at the shear peak was found out to the largest probed value 1 
on the CL 1701+6414 and CL 1641+4001 cases. 



16'. See Sect.|4]for more details 



be attributed to either tracking or focussing issues of the tele- 
scope. In most fields, no extreme outliers were present or could 
be easily identified. Only frames with average PSF ellipticity 
|(e)[ < 0.06 entered the coaddition. No clear distinction leaving 
a sufficient number of low-anisotropy frames was possible for 
the CL 1641+4001 and CL 1701+6414 fields. In these cases, all 
frames with |(e)| <0.10 were used for coaddition. 

We classify as galaxies all sources fainter than the brightest 
unsaturated point sources (r' AVT0 > r mi n ) an d more extended than 
the PSF (§>§* max ). Because even poorly resolved galaxies carry 



a lensing signal, we add sources r'. 



AUTO 



> r 



and & > ?? a 



with i?™ a as 0.95i?^ ax to the galaxy shape catalogue (cf. Fig. 4 
of Paper I). The parameters defining this catalogue for each field 
are tabulated in Table [2] together with its number density « ga i « 
"ksb/2. 

2.4. Background Selection 

Cluster weak lensing studies rely on carefully selected cata- 
logues of background galaxies, the carriers of the lensing sig- 
nal. While falling short of yielding a reliable photometric red- 
shift (photo-z) estimate for each individual galaxy, three-colour 
imaging makes possible a selection of foreground, cluster, and 
background sources based on their distribution in colour-colour- 
magnitude space (cf. |Medezinski et al.1|2010| Klein et al., in 
prep.). The method described below that we use for the three- 
band clusters is an improved version of the background selec- 
tion in Paper I, to which we refer for concepts and terminology: 
While considering all galaxies fainter than mf a i nt in the lensing 
catalogue, galaxies brighter than mbright are rejected. In the inter- 
mediate regime (weight < r' < fMf a intX we include galaxies whose 
g'-r' versus r'-i' colours are consistent with sitting in the back- 
ground of the cluster, based on the - similarly deep - Ilbert et al. 
(2009) photo-z catalogue (see Appendix [B] for more details). 

For the single-band clusters, the background selection sim- 
plifies to a magnitude cut, meaning that lensing catalogue con- 
sists of all galaxies r' > mf a j nt . For each of these sources, our 
KSB+ implementation yields a PSF-corrected ellipticity s = 
£1 +IS2, a noisy but, in principle, unbiased estimate of the re- 
duced gravitational shear g — g\ + ig2 (cf. |Schneider|2006| l. We 
choose the values of ;«f a i n t (and m^ght where applicable) such 



that the signal-to-noise ratio of the aperture mass estimator, or 
S -statistics ( Schneid er|1996| ) is optimised: 



S {0c\ #out) : 



V2Z;Stj!2/|0;~0clAW) 



OS 



^Lj Q){\0j-o c \ie oat ) 



(i) 



By s t j = Re[eexp(-2ii£>)], we denote the tangential ellipticity of 
the galaxy at position 8j, which with respect to the point C has 
a phase angle if. Equation ([T]) considers the noise from intrinsic 
source ellipticity, measured as o~ e = {e\+ e 2 ^ 1 ^ 2 ; while Qj(\0j- 
C \/O out ) is the |Schirmer et alT| ( |2007| l filter function with outer 
radius # out , maximising 5 for a cluster-like radial shear profile. 

We evaluate Eq. ([T]i on a regular grid with 15" mesh size. 
With the notable exception of CL 1701+6414 (Sect.|4T}, the sig- 
nal peaks are found close to the ROSAT-determined cluster centres 
and can be easily identified with our clusters. The adopted values 
of mbright, Wfaint, and 6 oul , yielding an optimal signal-to-noise ra- 
tio 5 max are summarised in Table|2] as well as the number density 
ni c in the lensing catalogu^] We refer to the grid cell in which 
Smax occurs as the cluster shear peak and discuss the significance 
of our cluster detections in Sect. 17.11 



2.5. Shear profile modelling 

Pursuing the approach adopted in Paper I, we model the tangen- 
tial ellipticity profile s t (8) of our clusters with the reduced shear 

~ [T996l |Wright & Brainerd] 
qi995[|1996||1997| 



profile g(6\ r2oo, cnfw) (jBartelmann 



2000) corresponding to the |Navarro et aL 
NFW) density profile. From the estimate of the radius r2oo - de- 
fined such that the density of the enclosed matter exceeds the 
critical density p c (zd) at the cluster redshift Zd by a factor of 
A = 200 - we infer the cluster mass M200 via 



A 47r t \ 3 

= AypcUd)r A 



(2) 



3 In addition to the photometric cuts, we restrict ourselves to high- 
quality sources defined by |£|<0.8, SExtractor detection significance 
v>4.5 and tr(P g )>0.1 for the KSB pre-seeing polarisability tensor. 
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Table 3. Additional parameters defining the "default" cluster models. Columns (1) to (5) tabulate the cluster redshift Zd, the average 
{{B}) and dispersion cr((J3)) of the distance ratio estimated from the CFHTLS Deep fields, the estimated fraction f& of residual 
foreground galaxies in the lensing catalogue, based on the same Ilbert et al. (2006 ) photo-z catalogue, and whether or not a dilution 
correction /i(6>) has been applied. Columns (6) and (7) present the celestial coordinates of the assumed cluster centre, with respect 
to which the fitting range r^ n <6<r max in columns (8) and (9) is defined. Columns (10) and (11) repeat the fitting range in terms of 
angular separation. Finally, we give the separation between rosat and lensing centres in column (12). 



Cluster 


Zd 


«yS» 


o-OJi)) 


/d 


Me) 


ttc,J20(X) 


<5c,J2000 


' min 


' max 




«max 


A6> 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


CL 0030+2618 


0.50 


0.348 r 


0.024 T 


0.152 T 


y 


OO^O"^^ 


+26° 18' 16" 


0.2 Mpc 


5.0 Mpc 


o:56 


14:04 


23" 


CL 0159+0030 


0.39 


0.447 


0.023 


0.087 


y 


01 h 59 m 18 s .2 


+00° 30' 09" 


0.2 Mpc 


5.0 Mpc 


0:65 


16:20 


79" § 


CL 0230+ 1836 


0.80 


0.168 


0.020 


0.321 


y 


02 h 30 m 26 s :6 


+18°36'22" 


0.2 Mpc 


5.0 Mpc 


o:46 


11:42 


20" 


CL 0809+28 11 


0.40 


0.437 


0.023 


0.105 


y 


08 h 09 m 41 s .O 


+28°11'58" 


0.2 Mpc 


5.0 Mpc 


0f64 


15:95 


178" § 


CL 1357+6232 


0.53 


0.324 


0.024 


0.198 




13 h 57 m 19i4 


+62°32'42" 


0.2 Mpc 


5.0 Mpc 


0:54 


13:62 


50" 


CL 1416+4446 


0.40 


0.437 


0.023 


0.136 




14 h 16 m 28 i :i 


+44°46'38" 


0.2 Mpc 


5.0 Mpc 


0f64 


15:95 


19" 


CL 1641+4001 


0.46 


0.381 


0.024 


0.127 




16 h 41 m 52 s .3 


+40°01'27" 


0.2 Mpc 


5.0 Mpc 


0:59 


14:70 


95" 


CL 1701+6414 


0.45 


0.381 


0.024 


0.134 




17 h 01 m 22 i;5 


+64°14'08" 


0.2 Mpc 


5.0 Mpc 


0:60 


14:88 


66" 



t We use an improved fit for the redshift distribution and corrected / d for CL 0030+26 18 with respect to Paper I. 
§ Lensing centre is within a large masked area, which probably diminished the accuracy. 



The best matching cluster mass profile parameters r2oo and cnfw 
minimise the merit function 



determine the fraction of residual foreground galaxies in the 
lensing catalogues (Table[3] cf. Sect. 6.1 1. 



X~ 



1=1 



\gi(0r, r 2 oo, c NFW )-s tJ (9i)\ 2 

^ 2 (ft)(i-te(ft;?-2oo,c NF w)l 2 ) 2 



(3) 



which we evaluate on a regular grid in r2oo and cnfw- By s t>( , 
we denote the tangential component of the scaled ellipticity 
Si = fofi(Oi)si for the z'-th galaxy, including a global shear cal- 
ibration factor /o = 1 .08 (see Paper I and |Hartlap et al.| [2009 ) 
as well as a separation-dependent correction f\(6) for the shear 
dilution by cluster members (detailed below). Accordingly , the 
error scales as &i(6i) = /o/i(#,)<x £ /V2 with cr E from Sect. 2.4 



The index i runs over all lensing catalogue galaxies with sepa- 
rations within the fitting range m i n <6< # max from the assumed 
rosat cluster centre, presented in Table[3] We choose separations 
# m j n and # max corresponding to distances of r m ; n = 0.2 Mpc and 
'"max = 5.0 Mpc at the respective cluster redshift. The denomina- 
tor of Eq. ([3| accounts for the dependence of the noise on 
itself ( [Schneider et al.|2000| >. 

Source redshift distributions The reduced shear g,(#,) exerted 
by a lens on the image of a background source further depends 
on the ratio of angular diameter distances between deflector and 
source and source and observer D s . For each of our fields, we 
estimate a catalogue-average (B) = (Ais/£ ) s)/=i...iV using the Ilbert 
et al. ( 2006 1 photo-z catalogue, drawn from the CFHTLS Deep 
fields with similar source number counts as a function of magni- 
tude r' as our MMT observations. Applying the same photomet- 
ric cuts as to the MMT data to the catalogues of reliable photo-z 
sources (cf. Paper I), we thus obtain proxy redshift distributions 
for our cluster observations. We repeat the fit of a van Waerbeke 
et al. (2001 1 redshift distribution and subsequent calculation of 
(B) as described in Paper I - but to an improved accuracy - for 
all combinations of MMT and CFHT Deep fields. As an input 
to Eq. Q, we use the mean {(fi))k=\..A (Table [3j measured for 
the Deep fields and consider its dispersion cr({B)) in the error 
analysis (Sect. [6*T| . 

We further employ the Ilbert et al. ( 2006 ) catalogue to test 
the efficacy of the background selection. Applying the respec- 
tive background selection to the Deep 1 photo-z catalogue, we 



Dilution by cluster members Although the selection of lens- 
ing catalogue galaxies is designed to include preferentially back- 
ground galaxies, we detect an increase in the fraction f rsc (0) of 
galaxies whose g' - i' colours are consistent with the red se- 
quences at Zd towards the centres of our three-band clusters. 
Tentative red sequence galaxies are defined using an interval in 
g' - i' empirically found in the g' - i' versus i' colour-magnitude 
diagram, around the expected colour of a Coleman et al. ( 1980) 



early-type galaxy calculated with the Bolzonella et al. (2000) 
photo-z code. To correct for the dilution effect of these likely 
unlensed sources in the shear catalogues, the corrective factor 
fi(6) = 1 + £(6>)/[£(#) + B] is introduced. The NFW surface mass 
profile £(#) and background term B are determined by a fit to 
/rsc($)- We apply this correction only to the three-band clusters 
for which the g'-i' information is available (see Table[5]). Because 
we have f\(8) measured for only four clusters, three of which 
suffer from large masks in the crucial central regions, we decide 
against using an averaged f\(0) for the single-band clusters at 
this stage of the survey. 



2.6. Surface mass maps 

While we use the tangential shear profile to determine cluster 
masses, we are interested as well in the (projected) mass dis- 
tributions of our clusters in order to distinguish possibly merg- 
ing systems of disturbed morphology from relaxed clusters. The 
non-local relation between shear and convergence K = £/£ cr j t can 



be inverted, as shown by Kaiser & Squires ( 19931. We perform 
mass reconstructions using the |Seitz & Schneider] ( |1996| |2001| ) 
finite-field inversion algorithm. Concerning the mass sheet de- 
generacy (cf. Schneider 2006 ), the mean k along the edge of the 
field-of-view is assumed to vanish. 

The dimensionless surface mass kock, with an arbitrary nor- 
malisation, is calculated on a regular grid. Because each cluster 
field has to be divided into an integral number of grid cells, the 
mesh size cannot be fixed to the same constant for all clusters, 
but varies slightly, with a mean of 40793 and a standard deviation 
of 07 36. For all clusters and grid points, the algorithm accounts 
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for lensing catalogue galaxies within a radius of S = 2'. The in- 
put shear field is smoothed with a truncated Gaussian filter of 
0.555 9 S full-width half-maximum, which drops to zero at 6> s . 



3. Results for normal clusters 

In this Section, we present the outcome of the WL modelling, 
by showing a comprehensive figure combining the lensing sig- 
nal maps, shear profile, and NFW modelling for each cluster. 
CL 1357+6232 (Fig. |2]i serves as our example; for more de- 
tails on the other clusters, we refer to Figs. |D.1| to Figs. |D.4| 
in Appendix [D] Two clusters, CL 1701+6414 (Fig. [3j and 
CL 1641+4001 (Fig. B), exhibit multiple shear peaks and shear 
profiles that are very flat but positive over a large radial range. 
The more involved modelling of these "special cases" - as op- 
posed to the "normal clusters" - is described in Sect. [4] 

In the upper panel of Fig. [2] we present the S -statistics 
(solid orange) and ^-contours (green dashed) for CL 1357+6232, 
overlaid on a cut-out of the Megacam r'-band image with ~20' 
side length. Masked areas can be identified from the red poly- 
gons (mostly squares). The rosat centre is given by a yellow, 
eight-pointed star symbol. A filled orange square denotes the 
shear peak grid cell (Sect. |2,4), while a star symbol with error 



bars shows the WL centre from bootstrapping (Sect. 3.2 1. 

The lower left panel of Fig. [2] shows the binned shear pro- 
file (s t (6)) as filled circles with error bars giving the disper- 
sion of the measured values. Open diamonds give the cross 
component (s x (9)) which is on average expected to be con- 
sistent with zero for cluster lenses. The red solid line denotes 
the best-fit NFW model. Finally, the lower right panel of Fig. [2] 
presents A;^ 2 (r2oo,CNFw) - X 2 ~ m ^ n (X 2 )- The minimum is indi- 
cated by filled circle; contour lines enclose the 99.73%, 95.4%, 
and 68.3% confidence regions, (i.e. Ay 2 =2.30, 6.17, and 1 1.30). 
An upward triangle marks the minimum of Ax 2 when restricting 
cnfw to its B ullock et al.| ( |200T| dashed line) value. 

3.1. Cluster detection and lensing morphology 

We successfully detect all observed 400d clusters using the S - 
statistics with at least 3.5<x significance and are able to derive 
a weak lensing mass estimate for each cluster. Table [2] sum- 
marises the maximum detection levels S and the optimal filter 
scales 9°^ v The most significant detection is CL 0030+26 18 at 
z = 0.50 with S =5.84 (Paper I); the formally least significant de- 
tection is CL 0230+ 1836 at z = 0.80 with S =3.64. The S =3.75 
measured for CL 1701+6414 has a contribution from the nearby 
cluster A 2246 at 9 w 270" separation (Sect. 4. 1 1, rendering it the 
least secure detection: For f? out = 220", we detect CL 1701+6414 
at the 2.5<x level. By detecting CL0230 +1836, we demonstrate 
the feasibility of Megacam WL studies out to the highest red- 
shifts accessible for current ground-based weak lensing. 

In general, we find a very good agreement between the sig- 
nal morphologies, of the S -statistics and mass reconstruction, 
i.e. we detect the same structures at comparable relative signal 
strength. This result reaffirms that our detections are not caused 
by artifacts in the (independent) analysis methods. 



3.2. WL cluster centres 

We define a "default" model for the NFW modelling of each 
cluster, determined by the parameters in Table [3] i.e. the cluster 
centre, fitting range, ((J3)), and dilution correction. We acknowl- 
edge that a careful and consistent treatment of cluster centres is 



important to prevent masses from being biased. In the default 
model, we use the lensing-independent rosat X-cluster centres. 
For comparison, we also consider cluster centres based on the 
S -map, which provide us with a high signal-to-noise shear pro- 
file. The shear peaks (most significant cell in the S -map) are 
thoroughly studied with respect to the background selection pa- 
rameters and their interpretation as significances (Sect. [77T| . 

The 5 -peak of CL 0159+0030 is located conspicuously 
close to the edge of an extended shear plateau which is likely 
caused by a large masked area [^] around a bright star (V = 8.3, 
Figs. [T] and D.l i. Similarly bright stars are present also close to 
CL 0230+ 1836 and CL 0809+28 11 (Figs. [112] and |D3J. In the 
latter case, where the S -peak lies within the masked area, we 



discuss the effect of masking in Sect. 7.5 

As noise can boost S in a grid cell compared to its 
neighbours, we perform a bootstrap resampling of the S -map 
(cf. Paper I) in two cases, CL 1357+6232 and CL 1416+4446. 
Averaging over 10 5 realisations, for which we draw N\ c galax- 
ies with repetitions from the lensing catalogue, we determine 
a lensing centre. We find the bootstrap lensing centres to be 
in good agreement with the shear peaks of CL 1357+6232 and 
CL 1416+4446, well within the standard deviation of the boot- 



strap samples. In Sect. 7.2 the implications of the choice of clus- 



ter centres for the mass estimates are discussed. 



3.3. Shear profiles and NFW modelling 

Five of our clusters can be classified as "normal", charac- 
terised by centrally increasing (s t )(9) profiles, in good agree- 
ment with the NFW models. As expected, their (s x )(9) pro- 
files are consistent with zero, with fluctuations that can be ex- 
plained by shape noise. The two other clusters, CL 1641+4001 
and CL 1701+6414 show a more complicated morphology in 
their S -maps (Sect. |4]l. 

Table [4] provides the cluster parameters resulting from the 
NFW modelling. Uncertainties in r2oo and cnfw are calculated 
from Ay 2 corresponding to a 68.3% confidence limit for one 
interesting parameter (Ay 2 = 1). Cluster masses M^^oo) are 



inferred via Eq. 



3.4. Mass-Concentration Relations 

Weak lensing hardly constrains cluster concentration parame- 
ters, because the dependence on cnfw is highest in the cluster 
centre where few lensed sources are observed. This is reflected 
also in our results, with huge uncertainties measured for Cnfw 
in several objects. Hence, we perform two additional measure- 
ments, in which we fix the value of cnfw- 

The first mass-concentration relation we assume is the one 



found by Bullock et al. (2001 B01) for simulated clusters 



cboi 



cboi,o / M, 



1 + z \ M m 



(4) 



with cboi.o = 9.0, a B oi = -0.13, and M, = 1.3 x 10 13 /r'M©. In 
their simulations, B01 observe a scatter of o-(logc v j r ) = 0.18 for 
a fixed M vir . 

For our purposes, we insert M^J (r2oo) for Mvir in Eq. 
Due to the weak dependence of cboi on M v ; r this results only in 
a very small underestimate of c B01 . For two of the total of eight 



4 Due to the filtering with large scales 9 oM , we measure a signal also 
in masked areas. Naturally, the correlation between neighbouring grid 
cells is even higher than in unmasked regions. 
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^2000 

13 h 58 m 30 s 13 h 58 m s 13 h 57 m 30 s 13 h 57 m s 13 h 56 m 30 s 13 h 56 m s 



62°40' 



62°24' 



V 0.05 



A 

«r o.oo 

v 



-0.05 



i — • • r . — » t v»< 

CL 1 357+ 6232 < 

— 




NFWfit 
binned <£„> 



centred on S-peak 



• free c 200 parameter 
a assuming c B01 
t assuming c B12 
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o 

Q. 



10 



15 



on 



10 



15 



'NFW 



Fig. 2. Lensing results for CL 1357+6232. Upper panel: Megacam r'-band image (cut-out of ~ 20' side length), overlaid with S- 
statistics (orange solid) and k (green dashed) contours. Contours start at S = 1.0 and k = 0.01, with increments of AS - 1.0 and 
Ak = 0.01, respectively. The eight-pointed star symbol marks the rosat cluster centre, while the filled square shows the shear peak 
grid cell. A star symbol with error bars denotes the lensing centre from bootstrapping. Regions inside red polygons are masked 
out of the analysis. Lower left panel: Profiles of the binned tangential «e t ), filled circles) and binned cross ((s x ), open diamonds) 
ellipticities. Error bars give the bin dispersions. The thick blue curve denotes the best-fit NFW model (Eq. Here, the S -peak 
was assumed as centre. Lower right panel: AY 2 (r2oo,CNFw) w ^ tn res P ect t0 its minimum (filled circle), shown as grey shading and 
contours indicating lcr, 2<x, and 3<x confidence. An upward (downward) triangle on a dashed (dash-dotted) curves mark the best-fit 
values assuming the B01 and B12 mass-C2oo-relations. A diamond marks cluster parameters assuming the S-peak as centre. 
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Table 4. Synopsis of cluster parameters and resulting weak lensing masses: Columns (1) to (3) give r2oo and cnfw fr om the default 
cluster models, and the corresponding mass, M^Li^m)- Fixing the concentration to the value cboi expected from the 



( j2001[ > mass-concentration relation, we obtain the best-fit radius t"2oo,boi and mass M^J OB01 (r2oo,BOi)> m Columns (4) to (6). Columns 
(/) to (9) contain the respective quantities calculated assuming the Bhattacharya et al. ( 201 l| l mass-concentration relation. 



Bullock et al. 



>"200,B01 t'BOl.NFW ^200,B01 ( r 200,B01 ) 



Cluster 



^200 
[Mpc] 
(1) 



cnfw 



(2) 



1 ^9+0.14 

1 -27+0.18 

-Q-J2 
1 «4+0.28 
1,J -0.32 

1 75 + - 23 
94+0.32 



1 7+ 13 ~ 
1 • ' -0.8 

>20 
5 5 +t ™ 

—3.3 

r - 3 

y >- 1 n 1 



[1O 14 M ] 
(3) 



7 25+2.W 

4 67+^ 
10.78!™ 
9.841^ 
2 92 +1 99 

-0-79 
9 o 4 +2.61 

-1.08 



[Mpc] 
(4) 



(5) 



1 aq+0.13 t q+2.0 
1-J -0.14 -1.3 
1.34+0-21 43+2.2 



[10 14 M o ] 
(6) 



5 50 +rar 

-1.50 
4 qq+2.44 

-2.03 
9 68+ 602 

OT+3.15 

2 26+' :4?i 
-1.3) 

1.69+o-B 

1 69+ 



^2003 12 
[Mpc] 
(7) 



t'Bi2.NF\v 



(8) 



200.B12 



(^200,812) 



1 45+OTT3 Q 4+I.7 

1 31+ - 22 3 6+ L9 

• J -0.27 J- -1.2 

1 c 1+0.27 9 n+1.5 

w -0.32 *" y -\B 

1 7 r+0.19 o 4+1.7 

-0-22 J ' H -1.2 

1 06+ 8 - 20 3 6+ 18 

" -0.24 -1.2 

n QQ+0.T6 t n+2.0 

-Q18 -1.3 

1 01+ 8 - 2 ' 3 7+ 1 ' 9 

l-0ll! 3.8+i| 



[10 14 M o ] 
(9) 



6 23+ L87 

"■■"-1.76 

4 10+ 246 

' -1,76 

10.23^=1 
9.24+ 346 
2 52 

1 71+0$ 

205^ 

-1-26 

2 03+ LT " 



CL 0030+26 18 
CL 0159+0030 
CL 0230+1836 
CL 0809+28 11 
CL 1357+6232 
CL 1416+4446 
CL 1641+4001 
CL 1701+6414+ 



1.62 



1 4Q+0-26 T Q+li 

1 69+°- 18 3 9+ 2 

102 +o.T8 4.2+ 2 - 2 

97+B:?I 4 g+ii 

0.86+II 4.S-J; 



t Fixing the radius and concentration of A 2246 to r Si 2oo = 0.90 Mpc and cnfw = 20 and hence using A^ 2 = 1 for the lcr error margins. 



clusters we analysed, cboi is very close to the cnfw obtained by 
lensing, while for others it differs strongly (see TableQ. 

Assuming the B01 mass-concentration relation, we 
apply a Gaussian prior p c (r2oo, Qoo) with standard de- 
viation cr(logC2oo) - 0.18 to the tabulated values of 
A^ 2 ( r 200> Cnfw) f° r eac h of our clusters, and marginalise over 
the C200 dimension. The radii r2oo,Boi and the correspond- 
ing masses M 2 ^ obo1 (/2oo,boi)> found from the minimum of 

Tjj Pcteoo, c 2 ooj)Ay 2 (r 2 oo, Qooj) are listed in Tabled 

We notice that the simulations from which the B01 rela- 
tion was measured assume <r 8 = 1.0 to fix the normalisation of 
the matter power spectrum. This value is inconsistent with more 
recent measurements of cosmological parameters (e.g. [Larson 
|et al.|[20TT] |Burenin & V ikhlinin 2012). Hence, we consider a 
second mass-concentration relation, based on a recent suite of 
simulations employing cr% = 0.8 as favoured by currents mod- 
els: |Bhattacharya et al.| ( [20TT] B 12) study Dark Matter haloes of 
massive clusters and find that the concentration parameter can 
be modelled with a single power law when expressed in terms of 
the peak height parameter v from linear collapse theory]^] Their 
simulated clusters are best represented by: 



C200,Bi2(v) = £>(z) a5 x5.9v- () - 35 



(5) 



with a variance of cr c = 0.33c2oo- We compute the growth fac- 
tor D(z) for a flat Universe with a cosmological constant. In 
complete analogy to the B01 relation, we compute cluster radii 
r200,Bi2 and masses M™ 1 OB12 (r 20 o,Bi2) for each cluster, given the 
C2oo-^20o-relation resulting from Eq. (|5j. The results are pre- 
sented in Table [4] 



4. Special Cases 

4.1. CL 1701+6414 

4.1 .1 . X-ray clusters and shear peaks 

A weak lensing analysis of CL 1701+6414 has to deal with shear 
by multiple structures. The strongest shear peak {S = 4.3<x) in 
Fig. [3] coincides with the most prominent cluster in the field 
amongst optical galaxies^] Abell 2246 (big four-pointed star 



We use the fitting formula for v(M, z) from Table 2 of Bhattacharya 



etal. 



1 20111. 



symbol in Fig.[3), 4'.2 to the west of CL 1701 +6414. With a red - 
shift of z = 0.225 ( [Vikhlinin et"aT1[T998l |Burenin et aLl|2007) , 
A 2246 is part of the 400d parent sample, but not of the distant 
cosmological sample. CL 1701+6414, for whose detection in the 
S -statistics the lensing catalogue was optimised, is detected at 
the 3.7cr level. The rosat catalogue of | Vikhlinin et al. (1998) 
lists two further clusters in the field, VMF 191 at z— 0.220 and 
VMF 192 at z=0.224 (small star symbols in Fig. [3j, which we 
identify with S -peaks of 2.9<x and 2.7cr, respectively. Another 
3. lcr peak lies close-by. A zone of positive shear signal extends 
over >20', from the north-east of VMF 192 to a 3.6<x shear peak 
south-west of A 2246, which does not correspond to a known 
cluster. Noticing the very similar redshifts of A 2246, VMF 191, 
and VMF 192, we likely are observing a physical filament at 
z - 0.22, through whose centre we see CL 1701+6414 in pro- 
jection. Luckily, a likely strong lensing arc, 10" to the west of 
the BCG of CL 1701+6414 (z = 0.44+0.01, |Reimers et al.|1997) 
gives direct evidence that CL 1701+6414 acts as a gravitational 
lens. We find no significant WL signal near the rosat source 
RXJ1702+6407 ( |Donahue et al.|2002| cf. Appendix [Tl6i. 



Due to the high concentration of galaxies, the region is masked. 



4.1.2. Two-cluster modelling of MMT data 

Plotting the binned tangential shear around the lensing centre 
(lower left panel of Fig. |5J, we find a flat profile whose av- 
erage (s t (6)) > is consistent with the extended shear signal 
in the S -map. The cross-component (s x (0)) is consistent with 
zero. Acknowledging the prominent signal related to A 2246, we 
model the shear of CL 1701+6414 and A 2246, simultaneously, 
using an NFW shear profile for each deflector. 

We assume both the shear g p of the primary and g s of the 
secondary component to be small. In this limit, the shear com- 
ponents originating from both lenses become additive: 

gadd.a (0) = gp,a (0', * p,200 . Cp,NFw) + gs,a(0\ r s,200 . Cs.NFw) , (6) 

with a - 1,2. Here, r p>2 oo, r s , 2 oo, c P;N fw, and c SiN fw are the 
radii and concentration parameters of the primary and secondary 
component. Note that g a dd,»(0) explicitly depends on the two- 
dimensional coordinate vector 6: the shear field of two clusters 
no longer has radial, but only axial symmetry. This is illustrated 
in the lower right panel of Fig. [3] showing the shear fit expected 
from the best-fit two-cluster model for the CL 1701+6414 lens- 
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Fig. 3. Shear signal in the CL 1701+6414 field and its best-fit model with two NFW components accounting for CL 1701 +6414 and 
A 2246. Upper plot: The layout follows Fig. [2] The rosat position of A 2246 is marked by a big four-pointed star symbol. Smaller 
star symbols denote positions of further X-ray clusters. Lower left plot: The layout follows Fig. [2] The solid blue and dashed red 
lines give the mean tangential and cross shear components, averaged in bins around the CL 1701+6414 shear peak, as expected 
from the two-cluster model. The separation of the two main clusters is indicated by a vertical dotted line. Lower right plot: The 
orientations and amplitudes of the shear, as expected from the best-fit two-cluster model, calculated on a regular grid. 
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Fig. 4. Simultaneous NFW mod- 
elling of CL 1701+6414 and 
A 2246. Each panel shows the 
dependencies between two of 
the four parameters, with the 
other two marginalised. Solid 
confidence contours (lcr, 2cr, 3<r) 
denote the default case, using the 
rosat centres; dashed contours 
denote models centred on the 
5 -peaks. The respective parame- 
ters minimising x\ Me indicated 
by a filled circle and a upward 
triangle. Sets of thin contours 
denote the confidence contours and 
parameters minimising x\ obtained 
from the analogous analysis of 
the CFHT lensing catalogue. The 
best-fit values for the CFHT data 
are marked by downward- and 
upward-pointing triangles for 
the ROSAT and 5 -peak centres, 
respectively. 



ing catalogue, evaluated on a regular grid. We consider a modi- 
fication of the merit function given by Eq. 



, v-i gadd,iVp,200, CnNFW, ^s,200; C s ,NFWl _ fi/ 
i=l 



(7) 



gadd,;(>"p,200, Cp.NFW, r s,200, Cs.NFw) 



The symbol x\ highlights the dependence on four parameters, 
the radii and concentrations of the two clusters. Note that x\ 
models the measured e, directly, without recourse to a definition 
of the tangential component. 

We assumed «/?» = 0.381 for CL 1701+6414 (Table [3} 
and ((/?)) = 0.640 for A 2246, calculated the same way as for 
the other clusters. The average tangential and cross-component 
of the shear expected in concentric annuli around the centre of 
CL 1701+6414 are presented in the lower left panel of Fig. [3] A 
vertical dotted line denotes the separation of CL 1701+6414 and 
A 2246. We find a good agreement to the measured shear and 
note that due to the lack of radial symmetry the dispersion of the 
model values in the annuli is of the same order as the measure- 
ment errors. Although the cross-component can be large at some 
points in the image plane, (g x ) cancels out nearly completely 
when averaging over the annuli. 

Figure |4] presents the confidence contours and parameters 
minimising Eq. |7]i for the default model (filled circle and solid 
contours). The panels of Fig. [4] show all combinations of two 
fit parameters, where we marginalised over the two remaining 



ones. Owing to the 4-dimensional parameter space, we tested 
a coarse grid of points to avoid excessive computing time. The 
picture emerges that r p 2oo and r Sj 2oo are relatively independent of 
each other (top right panel). Hence, the presence of the respec- 
tive other cluster does not seem to affect the accuracy with which 
we can determine the masses of the two clusters strongly. The 
data favour the smallest tested value, c Pj nfw = 0.5 for the concen- 
tration of CL 1701+6414, and the largest one, c Si nfw = 15.5, for 
A 2246. Using shear peak cluster centres (dashed contours and 
upward pointing triangle in Figure]?]), c p nfw is also very low, but 
the uncertainties are large. The poor constraint on c Si nfw might 
be partly due to the masking of the centre of A 2246 or shear 
contribution by the BCG. 



Given the absence of a strong covariance between the pa- 
rameters of A 2264 and CL 1701+6414, we fix the parameters 
of the foreground cluster to r S 2oo = 0.90 Mpc and c s NFW = 20 
and repeat the analysis with our usual, finer parameter grid. The 
best model is found for r™ = 0.94+° ;g Mpc c™ FW = 0.1+ 1 ,;!, 
confirming the results from Fig. [4] We note that we find a low 
Cp.nfw, although our model explicitly accounts for the extra 
shear by A 2246. Using the default model, we compute masses 
of 1.6+^ x i0 14 M o forCL 1701+6414 and l.l^xl0 14 M© for 
A 2246, based on r™ = 0.9 + 0. 1 Mpc. 
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Fig.5. Same as Pig. [2] but for CL 1641+4001. In the map, a triangle denotes the secondary shear peak, while a small star symbol 
marks the position of the |von der Linden et al. ( 2007) > cluster candidate. Note that no peak in the complex pattern of shear peaks 
correlates with its position. 
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4.1 .3. Comparison to CFHT data 

In addition, Fig. [4] shows confidence contours obtained from a 
WL analysis of CFHT observations of the CL 1701+6414 field 
(r'-band, w 7200 s), which we discuss in greater detail in Sect. [5] 
We repeated the two-cluster modelling using Eq. (|7]) follow- 
ing the same data reduction and shear measurement pipelines. 
Besides mf a i nt = 20.2 and the PSF-dependent galaxy selection, 
parameters are kept at their MMT values. 

The resulting cluster parameters minimising x\ (red down- 
ward triangles for rosat and diamonds for S -peak centres) and 
the corresponding thin confidence contours in Fig.|4]show agree- 
ment with the MMT data within the 2cr margins or better. With 
r™ n 00 = 0.70 ± 0.20 Mpc, and r™ = 0.85^ Mpc, relating to 

WL masses of M^™ T = 0.7+^xl0 14 M for CL 1701+6414 and 

M 200 HT = 09 ± 0.4x 1O 14 M for A 2246, we arrive at lower 
masses, especially for CL 1701+6414, but consistent within the 
uncertainties of the MMT data. Using the S -peaks as centres 
yields very similar results. 

Our CFHT data give more plausible best-fit concentration 
parameters of c™ FW = 7.5+> 8 for A 2246, and c^w = 4 - 5 -4S 
for CL 1701+6414, although the constraints are poor. We con- 
clude that a dual-NFW modelling is feasible, but more sensi- 
tive to the choice of cluster centres than a single NFW fit to f-200 
and c NFW . Adding more cluster components would even increase 
these interdependencies. However, the main point here is that the 
MMT and CFHT analyses agree. 

4.2. CL 1641+4001 

The S -statistics map of CL 1641+4001 exhibits several shear 
peaks which form a connected structure of >20' extent (Fig.[5]l. 
Located within a plateau of S >3cr significance, the rosat centre 
of CL 1641+4001 (big star symbol) is separated by 95" from the 
primary (S =4.12) shear peak and by 125" from the secondary 
(S - 3.95, orange triangle in Fig. [5]l shear peak. The BCG of 
CL 1641+4001 can be found between the rosat centre and pri- 
mary shear peak. 

The (s t (8)) profile (lower left panel of Fig. [5]) centred on 
the main shear peak profile is flat, with a positive average in all 
bins and the most significant positive signal at ~9' distance from 
the cluster centre. In the innermost two bins (#<3!33), (e x {8)) is 
of similar amplitude as the tangential component, but consistent 
with zero at the lcr level. Similar to CL 1701+6414, the mod- 
elling using Eq. (|3jl finds a very low cnfw = 0.1+qj, consistent 
with zero and reflecting the flat shear profile. 

The only cluster candidate in the literature besides 
CL 1641+4001 is SDSS-C4-DR3 3628 at z _ 0.032 identi fied 
in the SDSS Data Release 3, usi ng the |Miller et al.| (|2005|> al - 
gorithm, but published solely by von der Linden et al. { 2007) . 
We test a two-cluster model, introducing a second component 
at the redshift of SDSS-C4-DR3 3628, implying «/?» = 0.940. 
We choose the second-highest shear peak as the centre of the 
secondary component. The offset of ~ 3' to the coordinates of 
SDSS- C4-DR3 3628 (small star symbol in Fig. [5} is justified by 
the large mask at the latter position. The two-cluster fit yields a 
mass of order 10 14 Mq for both the primary and the secondary 
component. This estimate is in stark disagreement with the ab- 
sence of a massive, nearby cluster from our MMT image, which 
would have had to be missed by all but one cluster surveys. 

At the same coordinates as SDSS-C4-DR3 3628 and also 
at z = 0.032, NED lists CGCG 224-092, a galaxy pair, domi- 
nated by the bright elliptical UGC 10512. These two galaxies 



Table 5. Observation dates, final exposure times and seeing val- 
ues in the coadded CFHT/MegaCam data for CL 1701+6414. 



Filter 


Observation Dates 


T 

1 exp 


Seeing 


g' 


2006-03-07 


1601s 


0794 


r' 


2006-05-29 


7179 s 


a/66 


i' 


2006-04-26 


1922 s 


a/ 84 


z' 


2006-04-22, 2006-04-26 


1801s 


a' 82 



are what we see in the Megacam image^Jand also in the SDSS 
image of the area. Inspection of the respective Chandra im- 
age shows significant X-ray emission, whose extent of ~ 30" 
in diameter (~ 20kpc at z — 0.032) is consistent with being 
caused by a massive elliptical galaxy or small galaxy group. 
With w 1.7 x 10 4I ergs~' in the 2-10 keV range, its flux is 
high for a single galaxy, but the low temperature of a 0.6 keV 
(obtained by fitting an absorbed APEC model) speaks against a 
galaxy group. In conclusion, we deem it unlikely that the com- 
plex structure in the S -map of CL 1641+4001 bears a significant 
contribution from the z = 0.032 structure. 

We prefer the hypothesis that the shear is caused by a 
complex structure at the redshift of CL 1641+4001, although 
its X-ray morphology does not hint at a merger (Vikhlinin 
et al. 2009a I. Despite its shortcomings, we return to the sim- 
plest explanation for the time being and model CL 1641+4001 
by a single NFW component: We obtain a minimum of x 1 for 
+ 0.30™„„ n „ A - 0.1+0-3 These results, entail- 



min 1 Q(. 

'200 i-U -0.26 



Mpc and c 



min 

NFW 



-0.1 ■ 



ing a mass estimate of 2.3*?j x 10 14 Mq are illustrated by the 
filled circle and solid contours in the lower right panel of Fig. [5] 
Interestingly, choosing the secondary shear peak as a centre 
yields similar cluster parameters to those found by choosing the 
primary shear peak. This could hint at a major merger of simi- 
larly massive substructures, but more observations are needed to 
test this hypothesis. 



5. Verification with independent data 

5.1. CFHT Observations 

CL 1701+6414 is the only cluster we observed with 
MMT/Megacam for which deep, lensing-quality data ob- 
tained with another telescope exist. It has been observed in the 
g'r'i'z' filters (RL: G. Soucail, Run ID: 2006AF26) using the 
MegaPrime/MegaCam at the Canada-France-Hawaii Telescope 
(CFHT)0 Table |] lists the specifications of the CFHT data set. 
The CFHT data are processed with THELI in the same way as 
the MMT data, with a few CFHT-specific modifications to the 
code (cf. Erben et al. 2009), making use of the pre-processing 
available for archival CFHT data. Hence, the results are a suite 
of coadded and calibrated images in the g'r'i'z' passbands, 
centred on CL 1701+6414, and with a side length of ~ 1° each. 
From the g'r'i' images, we derived the pseudo-colour images of 
the centres of CL 1701+6414 and A 2246 in Fig.[T] 

We employ the CFHT data for two kinds of consistency 
checks with the MMT data: First, we run the lensing pipeline 
on the deep CFHT r' band image, applying the same shape re- 
covery technique to the same objects, but observed with dif- 
ferent instruments. The results of this comparison are detailed 



7 strongly overexposed and therefore masked 

8 For the sake of clarity, we use "MMT" and "CFHT" to distinguish 
the data sets. 
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Fig. 6. Sample density of the ratio e MMT / e CFHT f the complex 
ellipticities measured for the matched galaxies from the MMT 
and CFHT r'-band catalogues, respectively. The shaded con- 
tours correspond to the logarithmic densities of all galaxies 
from the MMT lensing catalogue which have a match in the 
CFHT catalogue. Solid contours give the density of galaxies 
detected with a signal/noise ratio of v > 15, the top 32.6 %. 
Note that the normalisation of the v > 15 galaxies is scaled up 
by 1/0.326 to obtain the same logarithmic contour levels. A 
Gaussian smoothing kernel of full-width half-maximum 0.075 
was applied to both contour maps. 



in Sect. 5.2 Second, making use of the CFHT imaging in four 
bands, we produced a BPZ (Benftez 2000) photometric redshift 
catalogue (Sect. 5.4 and Appendix |Q with the goal of test- 
ing the single-band (magnitude cut) background selection in the 
CL 1701+6414 MMT lensing catalogue. 



5.2. Comparative shape analysis 

In this Subsection, we compare shape measurements ob- 
tained in the same field (the one of CL 1701+6414) using the 
MMT/Megacam and CFHT/Megacam instruments (cf. Sect.|57T). 
Using the same parameter settings for our KSB pipeline, we 
extracted a KSB catalogue from the CFHT r'-band image. 
Subsequently, the CFHT and MMT catalogues were matched, 
using the associate and make_ssc tools available in THELI. 
With the smaller field-of-view of MMT/Megacam defining the 
location of possible matches, 68.2% of sources in the MMT KSB 
catalogue are matched to a CFHT detection. Larger masked ar- 
eas in the CFHT image - in particular due to reflections (so- 
called ghosts) around very bright stars - are the main cause im- 
peding a higher matching fraction. Inside the MMT area (mea- 
suring at a safe distance from its low-weight edges), we find 
85.5% of the CFHT sources to be detected by MMT. 

We note that objects in the matched catalogue have com- 
parable SExtractor signal-to-noise ratios v in both r'-band im- 
ages. Considering objects with vmmt > 15 - the top quartile of all 
objects in the catalogue of matches - for which selection effects 
should be negligible, we measure (vcfht/vmmt) = 0.832, with a 
dispersion of 0.057. These values show little dependence on the 
limiting value of v MMT , and confirm the visual impression that 
the r'-band images are of similar depth^ 

With these preparatory analyses in mind, we investigate the 
relation between the ellipticities observed with CFHT and MMT. 



9 While the good r cxp is similar for both data sets, the larger mirror 
area of MMT is probably offset by the better seeing in the CFHT image. 



Figure [^presents the ratio e MMr f s CFlirT f the complex elliptici- 
ties measured by KSB on the MMT and CFHT images{^]Shaded 
contours in Fig.|6]mark lines of equal density of the distribution 

MMT . r^Cl IT ' ' . . " . 



of e MMT /e CFHT , as measured from the sources passing the crite- 



ria for the MMT galaxy catalogue (cf. Sect. 2.4 1. Using a grid 



of mesh size as small as 0.01 for both the real and imaginary 



axes, we find the density distribution of g MM1 /g , - t ' m t scatter 
around its peak at unity. Note that the logarithmic scaling in 
Fig. [6] emphasises the wings of the distribution. When repeat- 
ing the analysis restricted to galaxies detected with vmmt > 15 
- the top 32.6 % of the matched sources contained in the MMT 
galaxy catalogue - the peak at g MMT = e CFHT persists, while the 
scatter is slightly reduced (solid contours in Fig. [6). This can be 
seen comparing the two outermost solid contours to the shaded 
contours, indicating the same levels of number density. 

This means, any systematic bias between shear measure- 
ments obtained with MMT and CFHT is smaller than a few per- 
cent. We expect a small bias, below the sensitivity of our mea- 
surement, to be present because of the dependence of the shear 
calibration factor fo on magnitude and half-light radius § (cf. 



Appendix C of Hartlap et al. 2009 1. In addition to our results 



from Paper I, the consistent galaxy ellipticities measured with 
MMT and the well-established CFHT/Megacam mark further ev- 
idence that MMT/Megacam is well-suited for measuring weak 
gravitational lensing signals. 

5.3. S -statistics from CFHT and MMT 

Figure [7] provides a qualitative comparison of the S -maps for 
CL 1701+6414 obtained with both CFHT and MMT. Its upper 
two panels show the independent shear catalogues drawn from 
the r' images of both instruments, at the respective optimal val- 



10 As the complex ellipticity is the relevant observable, we prefer con- 
sidering the components of the ratio £ MMT / e CFH T over me rat j os f or the 
individual components as measured with the two instruments. 



14 



Holger Israel et al.: The 400d weak lensing survey II 



MMT 




CFHT -> 



MMT matched 



10' 




CFHT matched 




Fig. 7. S -statistics in the CL 1701+6414 field drawn from the MMT (top left), CFHT (top right), and matched sources catalogues 
(bottom panels). The linear colour scale, contours indicating levels of S = 1 to S = 4, 6 out = 14'.5, and cross-hairs at the position of 
CL 1701+6414 are the same in all panels. Thick black squares outline the MMT field-of-view. 



ues m£™T = 21.9 and mf™ J = 20.2 for 6 0Ut = 14:5. The dis- 
tribution of the S -signal in the overlapping region inside the 
MMT field-of-view (black square in Fig. [7J) is astonishingly sim- 
ilar: Not only do we find the tentative filament from the north- 
east of VMF 192 to the south-west of A 2246 (compare Fig. [3] 
and the black lines in Fig. [7] indicating the aj2ooo an d feooo of 
CL 1701+6414). Moreover, also the regions of high S at the 
eastern and north-western edges of the MMT field-of-view cor- 
respond to peaks in the CFHT S -map. Whereas the detection 
significance at the peak closest to the position of CL 1701+6414 
is smaller for CFHT (S = 2.89 compared to S — 3.75), it is also 
more prominent in the sense of a deeper "valley" separating it 
from the dominant A 2246 peak (S =4.30 in both the MMT and 
CFHT maps). 



The second-most significant (4.08cr) shear peak in the 
CFHT S-map is at a J2 ooo = 17 h 01 m 57 s , foooo = +63°51', out- 
side the southern edge of the MMT field-of-view, with no known 
cluster but several brighter (/ < 20) galaxies in the vicinity. 

Can the subtle differences between the MMT and CFHT 
S -maps be attributed to shape noise or rather to selection of 
galaxies at the faint end? We investigate that by considering the 
matched-sources catalogue from Sect. |5.2| and apply to it the 
combined selection criteria for the MMT and CFHT lensing cat- 
alogues (e.g. both le^l < 0.8 and |e CFHT | < 0.8). The resulting 
S -maps derived from the MMT and CFHT ellipticities of the 
exact same sources are displayed in the lower left and lower 
right panels of Fig. [7] Naturally, all matched sources are lo- 
cated within the MMT field-of-view. Qualitatively, the matched 
S -maps again show the same structure, although they do not ap- 
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Fig. 8. Average and difference aperture mass At* (Eq. |8j, mea- 
sured in the matched MMT-CFHT catalogue. The grey-scale and 
white contours give AL ', the thicker blue contours show AL . 
The spacing for both contours is in multiples of 0.015, start- 
ing at 0. Black circles mark the positions of CL 1701+6414 and 
A 2246. 

pear to be much more similar than the S -maps drawn from the 
individual catalogues. This indicates that galaxy selection plays 
a relevant role. We note that based on the CFHT shapes of the 
matched galaxies, CL 1701+6414 is the most significant detec- 
tion with S =3.46, while the A 2246 peak is suppressed by « lcr 
compared to the pure CFHT map. 

Quantitatively, the Pearson correlation coefficient of g = 
0.912 between the matched-sources S maps substantiates the vi- 
sual impression of a high correlation. When the faintest « 20 % 
of galaxies are removed from the matched catalogue, consider- 
ing only galaxies brighter than r' cut = 24.85 in both the MMT and 
CFHT images, this value increases to g = 0.926. Removing the 
faintest * 40 % of galaxies by imposing r' cut = 24.45, it further 
rises to £ = 0.938. 

The detection of the same shear peaks reassures us that 
the multi -peaked S -distribution analysed in detail in Sect. |4.1| 
traces an actual shear signal and removes any doubts that the 
S -filament across the MMT field-of-view could be merely an 
instrument-dependent artefact, e.g. residuals of improper PSF 
anisotropy correction. 

As a final test to the hypothesis that we see the same shear 
signal measured in the ellipticities from both instruments, we 
consider the average and difference aperture mass of the matched 
sources: 

M ± =n6l uX N{0 c T l Yj <jQj(\ej-0 c \/0 oul ) (8) 



with sfj the tangential component of 



(£ MMT ±g CFHT )/2 



(9) 



While AL (blue contours) retrieves the signal of both clusters 
(black circles), exhibiting the expected great similarity to the 
matched-sources S -maps in Fig. [7] AL (grey-scale and white 



contours in Fig 
not obviously re 



8]) has a much smaller amplitude. Its pattern is 
ated to the one seen in AL: Although the main 
clusters reside in a region of enhanced At - , they do not corre- 
spond to peaks in At - . The absence of the AL -peaks in AL is 
consistent with the absence of a noticeable shear calibration bias 



between MMT and CFHT (cf. Sect. 5.2 1. A possible explanation 
for the stripe-like pattern in Ai~ are differences in the spatially 
varying anisotropy correction. We conclude these effects to be 
small and no impediment to direct comparisons of MMT and 
CFHT WL measurements, which we conclude to be consistent. 



The SW peak A peculiar feature in the MMT S -map of 
CL 1701+6414 is the 3.6cr shear peak at a mm = 17 h 00 m 05 s , 
£j2000 = +64°11'00), south-west of A 2246 (cf. Fig. |3). This 
peak does not correspond to an evident overdensity of galax- 
ies in the MMT and CFHT r'-band images, nor to an extended 
emission in the Chandra X-ray images. The pure CFHT S -map 
does not show a counterpart to the S W peak detected with MMT, 
although the S -contours of A 2246 are extended towards its di- 
rection. Interestingly, in the matched-sources S -maps, we de- 
tect a 2.8cr peak from the MMT ellipticities and a 2.4<x peak 
from the CFHT ellipticities. We checked that the SW peak in 
the MMT S -map does not arise from a chance alignment of a 
few galaxies with extraordinary high e t resulting from stochas- 
tic shape noise. Considering the observations in the CFHT and 
matched-catalogue S -maps, it seems likelier that we observe a 
true "shear peak" arising from the superposed light deflections 
of line-of-sight structure in the complex and dense environment 
of the CL 1701+6414/A 2246 field. 



5.4. Photometric redshift results 



Using photo-zs based on the available g'r'i'z' observations 
(Table [5} is challenging because of the small spectral cover- 
age and shallowness of the data. Nevertheless, a comparison 
with sources for which SDSS spectroscopic redshifts are known, 
revealed a coarse redshift sorting to be possible, with typi- 
cal errors of <x( z p h) x 0.25 for the relevant z < 0.5 redshift 



range (Appendix C.l i. Matching the photo-z catalogue with the 
CL 1701+6414 MMT data set, we can identify most sources in 
the galaxy catalogue with a photo-z galaxy, albeit with low qual- 



ity for most sources (Appendix C.2 1 



and the index j running over all N(0 C ) galaxies within a distance 
#out from C . The outcome of this experiment is shown in Fig. [8] 



Drawing the S -statistics from this catalogue, whereby 
galaxies are sorted based on their CFHT photo-zs, we retrieve 
shear peaks similar to Fig.|3]for the background catalogue, while 
CL 1701+6414 does not show up as a shear peak in the fore- 
ground catalogue ( Appendix |C.3| l. From this, we draw two con- 
clusions: First, we likely see an indication of shear emanating 
from more than one lens plane, namely CL 1701+6414 on the 
one hand and A 2246 and associated structures on the other hand. 
Second, owing to the poor quality of the four-band photo-z, a 
desirable calibration of single-band shear catalogues lies beyond 
the grasp of this data set. 
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Table 6. Weak lensing masses resulting from our analysis. Given are the WL masses M^j 812^200*12)' assumm g the B12 mass- 
concentration relation (cf. Table |4|, their lower and upper statistical (o"~ tat and cr^ tat ), systematic (tr sys and cr+ ys ), and total error 
margins (cr~ t and <x+ t ). In addition, the corresponding relative errors are presented. All masses are given in units of 1O 14 M0. 



Cluster 


MwlV^OO.wl) 


°Ttat 


<tat 




< ys 


o-tot 


<„, 


""stal 
M wl 


M wl 


°" sys 
M wl 


< 


M wl 


M wl 


CL 0030+26 18 


6.23 


1.76 


1.87 


2.07 


1.97 


2.72 


2.72 


28% 


30% 


33% 


32% 


44% 


44% 


CL 0159+0030 


4.10 


2.07 


2.46 


1.25 


1.31 


2.42 


2.78 


50% 


60% 


30% 


32% 


59% 


68% 


CL 0230+ 1836 


10.23 


5.18 


6.59 


4.68 


3.59 


6.98 


7.51 


51% 


64% 


46% 


35% 


68% 


73% 


CL 0809+28 11 


9.24 


3.14 


3.46 


2.59 


2.62 


4.07 


4.34 


34% 


37% 


28% 


28% 


44% 


47% 


CL 1357+6232 


2.52 


1.34 


1.74 


0.99 


0.89 


1.67 


1.96 


53% 


69% 


39% 


35% 


66% 


78% 


CL 1416+4446 


1.71 


0.79 


0.95 


0.60 


0.58 


1.00 


1.12 


46% 


56% 


35% 


34% 


58% 


65% 


CL 1641+4001 


2.05 


1.26 


1.51 


0.78 


0.73 


1.48 


1.68 


61% 


74% 


38% 


36% 


72% 


82% 


CL 1701+6414 


2.03 


0.98 


1.17 


0.78 


0.70 


1.25 


1.37 


48% 


58% 


38% 


35% 


61% 


67% 



Table 7. Components of the statistical error, assuming the B12 mass-concentration relation. We list all components entering (Eq. 10 1: 
The uncertainties cr^. due to shear calibration, and cr± eom from (D(j s /D s ), the projectional uncertainty cr* . due to cluster triaxiality, 

and cr^ ss from the projection of unrelated LSS. All errors are given in units of 10 14 Mq; the numbers in parentheses present the 
relative uncertainties. The relative statistical uncertainties are relatively independent of the choice of the concentration parameter. 



Cluster 


°~cali ""cali 




<r-. <r+ 

pro, pro| 


*^"lss *^"lss 


CL 0030+26 18 
CL 0159+0030 
CL 0230+1836 
CL 0809+28 11 
CL 1357+6232 
CL 1416+4446 
CL 1641+4001 
CL 1701+6414 


1.09(18%) 0.26(4%) 
0.38(9%) 0.14(3%) 
3.75 (37%) 0.58 (6%) 
1.19(13%) 0.39(4%) 
0.58 (23%) 0.14(6%) 
0.28 (16%) 0.07 (4%) 
0.41(20%) 0.12(6%) 
0.42(21%) 0.12(6%) 


0.65 (10%) 0.73 (12%) 
0.27 (7%) 0.29 (7%) 
1.56(16%) 2.36(23%) 
0.74 (8%) 0.80 (9%) 
0.26 (10%) 0.33 (13%) 
0.11(7%) 0.12(7%) 
0.24 (12%) 0.27 (13%) 
0.21 (10%) 0.21 (10%) 


0.62(10%) 1.00(16%) 
0.41 (10%) 0.66 (16%) 
1.02(10%) 1.64(16%) 
0.92(10%) 1.48(16%) 
0.25 (10%) 0.40 (16%) 
0.17(10%) 0.27(16%) 
0.21 (10%) 0.33 (16%) 
0.20 (10%) 0.32 (16%) 


1.52 (24%) 1.62(24%) 
1.08 (26%) 1.08 (26%) 
2.08 (20%) 2.08 (20%) 
1.97(21%) 1.97(21%) 
0.71 (28%) 0.71 (28%) 
0.44 (29%) 0.44 (29%) 
0.58 (29%) 0.58 (29%) 
0.58 (29%) 0.58 (29%) 



6. Accuracy of the mass estimates 

6.1. Error Analysis 

The error analysis of the seven clusters analysed in Sect. [3] fol- 
lows the method described in Paper I, i.e. we apply 

°it = °itat + °~sys = °"stat + ^LSS + °"pmj + °"geom + ^Lli ( 1 °) 

to calculate the total uncertainty in mass for each cluster. We will 
now discuss how we obtain the different terms in Eq. (fTOb. The 
statistical error cr stat is inferred from the tabulated for the 
cluster on the grid in r2oo and c NFW : Taking A^ 2 = 1, we find the 
upper and lower limits of r2oo an d then applying Eq. Q. Table[6] 
compares the masses of our eight clusters and their errors. 

The components cr ca i; and <r geom , accounting for the uncer- 
tainties in the shear calibration factor / and the redshift distribu- 
tion of the source galaxies are likewise determined from the anal- 
ysis of the parameter grid. Assuming the redshift distribution to 
be well modelled by the fits to the CFHTLS Deep 1 photo-z cat- 
alogue, we vary ((J3)) by the uncertainties tabulated in Table [3] 
As expected, cr geom increases with redshift because of the higher 
relative uncertainty in ((J3)). 

6.2. Redshift Distribution 

Comparing the source number counts in the CFHTLS Deep 1 
field with our MMT data, we find very good matches to the r'- 
band source counts in the CL 0030+2618 and CL 1641+4001 
fields, our observations with the deepest limiting magnitudes 
and a high density «ksb <; 40 arcmin~ 2 in the KSB catalogues 
(cf. Tables [T] and [2]). The other cluster catalogues exhibit a 



completeness limit (peak in the source count histograms) at 
slightly brighter r' -magnitudes, but follow the Deep 1 closer than 
the corresponding, alternative r + -band source counts from the 
COSMOS photo-z catalogue pibert et al.|2009[ l. In order to test 
for a possible bias in {/}) for the shallower cluster fields, we re- 
peat the fit to the redshift distributions from the four Deep fields 
with the following modification: Introducing a magnitude cut, 
we remove all galaxies with r' > r' m!lx from the CFHTLS cata- 
logues. Virtually independent of Zd, we find the (fi) for the cases 
with and without magnitude cut to agree within mutual error bars 
for rj nax >25.2, meaning that the variation within the Deep fields 
has the same amplitude as the effect of removing the faintest 
sources. In our shallowest field, CL 0230+1836, we measure a 
limiting magnitude of rf im = 25.1, with 15% of galaxies in the 
galaxy shape catalogue at r' > 25.2. We thus conclude that no 
significant bias in (fi) is introduced by using the/K// [libert et al." 
(2006 ) catalogue as a redshift distribution proxy and the disper- 
sion among the four fields as its uncertainty. 

6.3. Dilution by cluster members 

As in the case of CL 0030+26 18, we not only consider the un- 
certainty of ±0.05 we estimate for /o, but also take into account 
the dilution by remaining foreground galaxies in the shear cal- 
ibration error. Once again using the CFHTLS Deep 1 photo- 
Z catalogue as a proxy, we determine the fraction of galaxies 
at z p h < z c i after applying the respective background selection. 
We measure this fraction /d to increase with z: it varies from 
8.7% for CL0159+0030 to 32.1% for CL 0230+1836 (Table [3). 
As can be seen for CL 0809+2811 and CL 1416+4446 at the 
same redshift z = 0.40, the background selection based on three 
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bands results in a lower = 10.5% that the mere magnitude cut 
(/d = 13.6%) for only one band. Adding the two components of 
the error in quadrature, the lower limit we consider for fa ranges 
from 0.97 for CL0159+0030 to 0.73 for CL 0230+1836. 

6.4. Uncorrected Large Scale Structure 

To calculate the error <x L ss induced by LSS, we need to extrap- 
olate the findings of Hoekstra (2003), covering only the cases of 
5h~ 1 , 1 0/i~ 1 , and 20/T 1 xl0 14 M Q to lower masses. (Note that our 
M wl estimate for CL 0030+26 18 is very close to the f irst case.) 
The respective error contributions read from Fig. 6 of Hoekstra 
( |2003| l are ~ 1.2/T 1 , ~ 1.7/T 1 , and ~ 2.7/T 1 x 1O 14 M . By as- 
suming that the relative LSS error <x LSS /M wl increases linearly 
towards smaller masses, we arrive at the following relation: 

a- L ss/(10 I4 M o ) = fl M 14 + W< , (11) 



where a = 0.22/i , 
understand Eq. (11 



b = -0.01, and M u = M wl /(10 14 M©). We 
as an order-of-magnitude estimate for the 
LSS error and stress that simulated WL measurements are re- 
quired to provide a better understanding of this important source 
of uncertainty. In particular, we expect a larger cr LSS for higher 
Zd clusters, for which the existence of intervening massive struc- 
ture is more likely. We notice that the results of Hoekstra ( 2003| > 
are obtained at z = 0.3, more nearby than our clusters. 

In the special case of CL 1701+6414, the obvious LSS stz~ 
0.22 was taken into account by explicit modelling, in addition to 
what is described here. We estimate the error associated with the 
covariance of the parameters describing the two clusters in our 
simultaneous model for CL 1701+6414 and A 2264 (cf.|4~T| by 
varying the otherwise fixed r^oo of A 2264. The additional error 
from considering r s 200 = 0.8 Mpc and r Sj 2oo = l.OMpc, according 
to the uncertainties from our four-parameter fit, give a negligible 
contribution to the statistical error. Nevertheless, we caution that 
an additional uncertainty likely arises from our model choice. 

6.5. Triaxiality Projection Bias 



Applying the Kasun & Evrard (2005) fitting formula for the 
largest-to-smallest axis ratio of a triaxial halo as a function of 
mass to all our eight clusters, we arrive at expectation values of 
0.60 < rj < 0.64 for the largest-to-smallest axis ratio. Hence, con- 
sidering the triaxiality biases of |Corless & K ing (2007 ), we use 
cr + . = 0. 16 M wl for the error M wl induced by overestimation and 

proj ■> 

cr~ . = 0.10M wl for the one induced by underestimation caused 
by the projection of triaxial halos. 



7. Discussion 

The statistical, systematic, and total errors for all eight clusters 
are summarised in Table [6] both as absolute masses and as rel- 
ative errors. Table [7] provides the details on the composition of 
the systematic error for the eight clusters. We note that for all 
our clusters, in particularly the ones with small WL masses, the 
statistical uncertainties are the largest component in the total er- 
ror (the second largest usually being the projection of unrelated 
LSS). The relative statistical errors range between w 30% and 
~60%. The reason for this can be twofold: First, the large statis- 
tical uncertainties per se are caused by the small signal-to-noise 
in the lensing signals and thus a consequence of the low net expo- 
sure times in the lensing-band images, once we removed frames 
with high PSF anisotropy (Table [T}. Second, our account of the 



systematics might underestimate or neglect contributions to the 
systematic error. 

For instance, uncertainties in the determination of the cen- 



tres and the radial fitting ranges are not considered in Eq 
which we discuss in Sect. 



7.2 



fro), 

va- 



After checking the statistica 
lidity of our cluster detection in Sect. |7.1 we evaluate several ef- 
fects influencing the accuracy with which we can measure weak 
lensing masses. A reliable quantification of all of these uncer- 
tainties is beyond the scope of this article but very desirable with 
respect to the constraints on cosmology at which the next gener- 
ation of cluster weak lensing projects is aiming. 



7. 1 . Significance of Cluster Detections 

The S -statistics is known to produce spurious shear peaks even 
at high significance levels, although as simulations show, false 
detections above the ~4cr level are rare (e.g., Hetterscheidt et al. 



|2005[ [Dietrich et al.|2007[ ). Still, in principle, there is a nonzero, 
but small chance for one or the other of our detections to be 
false. Spurious detections are more sensitive against changes in 
the lensing catalogue or 6 out . Our tests with different photomet- 
ric cuts and values for 6> out found our cluster shear peaks to be 
robust. Another reaffirmation is the persistence of signals when 
bootstrapping the lensing catalogue, which we performed for 
CL 1357+6232, and CL 1416+4446. 

In order to test the interpretation of S -values as signifi- 
cances, we conducted the following test: For each galaxy in the 
catalogue, we add to the phase <p of the complex ellipticity esti- 
mator e = |e| exp (2i<p) an additional term ^ m( j drawn randomly 
from a uniform distribution in the interval < tp rn d < n. This 
procedure should completely remove the lensing signal from the 
data such that the resulting value of S be normally distributed 
around zero, with a standard deviation cr = 1 . 

We produced 10 6 realisations of such a randomised cata- 
logue for each cluster and find the S -distributions for all eight 
cluster detection to be well represented by a Gaussian distribu- 
tion. In all cases, the absolute of the mean value p of the fitted 
Gaussian is \p\ < 0.002, and of the same order of magnitude as 
the uncertainty in p derived from the fit. We do not find a bias 
to either positive or negative S . For six of the eight clusters, we 
find for the standard deviations cr of the fitted Gaussians values 
of I l-o-| < 0.01, with 0- = 0.962 and CL 1416+4446 and cr = 0.978 
for CL1641+4001, respectively, the largest measured deviations 
from the expected cr = 1 . 

For only one cluster, we find one \S\ > 5 event among 
the 10 6 realisations, consistent with the expectation of one such 
event in 1.7 x 10 6 realisations of the expected Gaussian distri- 
bution G. Therefore, we conclude that this randomisation test 
does not find indications for an overestimation of the signifi- 
cance of our cluster detections, as inferred from the S -statistics. 
On the contrary, the small standard deviations measured from 
the fits to the CL 1416+4446 and CL1641+4001 correspond to 
very slightly underestimated significances of these two cluster 
detections. 



7.2. The Role of Cluster Centres 



Dietrich et al. (2011! recently demonstrated that the S -peak 



gives a robust determination of the cluster centre, showing little 
susceptibility to projected large scale structure (LSS). However, 
using S -peaks as cluster centres for WL mass estimates is likely 
to result in a systematic overestimation of cluster masses, be- 
cause we pick such centres that produce the highest masses. 
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Fig. 9. Dependence of weak lensing cluster masses on the 
choice of centre and mass-concentration-relation. Upward and 
downward triangles show the masses using the B01 and B12 
mass-concentration-relation relative to the case of a free C2oo- 
The square and small triangles present models centred on the 
S-peaks for the free-C2oo, B01, and B12 cases, respectively. 
Clusters are ordered by increasing (free) C200 m the S -centred 
model. Grey lines denote the statistical error bars in the refer- 
ence model (RosAT-centred, free C2oo)- 



In addition to our default Rosat centres, we partnered all clus- 
ter models with a model choosing the S -peak as cluster centre. 
While we plan to quantify the centring bias once a more com- 
plete sample of our clusters is available, we can show overall 
trends based on the eight clusters discussed here. 

The separations between the shear peaks and rosat cen- 
tres are < 3' in all cases and < 1' for four of the eight clusters 
(Table [3]). This coincidence of X-ray and lensing centres adds 
further significance to the S -detections. Two out of the four re- 
maining clusters, CL0159+0030 and CL 0809+2811, have their 
S -peaks within larger masked areas, reducing the accuracy with 
which the centres can be determined. The complicated shear 
fields in the vicinities of CL 1701+6414 and CL 1641+4001, 
with separations > 1' between lensing and rosat centres are dis- 
cussed in Sects. l4~T1 and l4. 21 We note that the rosat cluster cen- 
tres themselves are accurate to ~ 10". 

First, we notice that the S -centred masses are biased high 
for six of the eight clusters (squares in Fig. [9]). The median is 
a 9 % higher mass using shear peak centres (^200™) tnan 
Rosat centres (Mf ™). We measure (M^ en /Mf m cen ) = 1.25, 
with a 0.36 standard deviation. This indicates the expected pres- 
ence of a centring bias, although a small one compared to its 
scatter. 

The most extreme differences are measured for the low- 
concentration clusters CL 1641+4001 and CL 1701+6414. The 
sorting in Fig.|9]by increasing C200 m the 5 -centred model seems 
to suggest a trend of a larger M2 00 cen /-^2(X) Cen in clusters of low 
concentration. We do not see a similar trend when sorting clus- 
ters by Af9, the separation between the two types of centres (cf. 
Table [3]), but caution that we deal with small number statistics. 
When applying the B01 (B12) mass-concentration relations (big 
and small upward and downward triangles in Fig. [9j, we arrive 



at similar results, with higher median biases of 18% (20%), but 
still large scatter. 

7.3. Low concentration clusters 

As Fig. [9] shows, the relative differences between masses mea- 
sured with the six different models for centre and concentra- 
tion parameter are most pronounced for CL 1641+4001 and 
CL 1701+6414, the two clusters with the lowest intrinsic con- 
centrations CjQy 6 ". 

First we remark, that low values for C200 are neither unex- 
pected nor unheard of in the weak lensing literature (cf., e.g., 
|Oguri et al.|20T0l [Planck Collaboration et al.|2012) . Small con- 
centration parameters as in CL 0809+28 1 1 or CL 1 357+6232 are 
not unique to our analysis method, especially when using a clus- 
ter centre significantly offset from the shear peak or in fields 
affected by masks in the data. Our analysis highlights the impor- 
tance of such effects. 



We note that, using an elliptical NFW cluster model, Oguri 
et al. ( 2010) > found four out of the 25 X-ray selected clusters they 
analysed a best-fit c v ; r < 1 . These four extreme cases are among 



the clusters Oguri et al. (2010 1 exclude from further analysis for 
centring problems or obvious misfit of the assumed model. The 
authors do not, however, report similar problems for three further 
clusters with 1 < c wlI < 2. 

The more extreme cases of CL 1641+4001 and 
CL 1701+6414 point to the limitations of describing a cluster 
shear field using a single, spherically symmetric function, even 
if obvious line-of-sight structure like A 2264 is taken into 
account. Nevertheless, tests like the comparison of Ar-maps for 
CL 1701+6414 from MMT and CFHT (Sect. [53) demonstrate 
that the measured low concentration parameters reflect the 
reality of the mass distribution and are not artifacts of the 
analysis. Concerning cosmological applications, the rejection 
of such clusters will likely result in a biased mass, which we 
prefer to avoid. We notice that, making different choices in the 
data analysis, High et al] ( |2012[ ) likewise succeed in measuring 
masses of clusters showing a flat shear profile. 

7.4. Influence of mass-concentration-relations 

Comparing the masses obtained with the three different choices 
for C200 (free parameter, B01 and B12 mass-concentration re- 
lations; Fig. [9] and Table [4), we observe interesting trends. 
Intrinsically, the B12 relation yields lower concentrations for the 
same cluster than B01. This relates to the result that the mass as- 
suming B01 exceeds the B12 mass, except for CL 0159+0030. 
Despite the overall agreement of the masses within their errors, 
the spread between the two masses seems to decrease with C2oo> 
independent of the choice of centre. 

On average, using a mass-concentration relation results in a 
lower mass estimate compared to a free C200, with ~ lcr signifi- 
cance for the B01 relation «M 2 oo,boi /-^200,def) = 0-85 ± 0. 16) and 
less so for the B12 relation «M 2 oo,Bi2/A^200,def) = 0.95 + 0.13). 

For the S -peak-centred masses, we obtain a similar pic- 
ture. The B01 masses increase relative to the free-C2oo masses, 
following this trend monotonically despite the statistical uncer- 
tainties. In Fig. [9] this can be seen looking at the squares and 
small upward triangles. The correlation can be explained by the 
shape found in the confidence contours of most of the clusters 
for which C200 can be constrained (e.g. CL 1357+6232; Fig. [2): 
Relatively high values for C200 oni y agree with the data for r2oo 
smaller than the best-fit value. The underlying reason is the ab- 
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Table 8. Three bright stars exacerbating the analysis of 
the CL0159+0030, CL 0230+1836, and CL 0809+2811 fields, 
identified by their BD and HD designations. We cite SIMBAD 
(http://simbad.u-strasbg.fr/simbad/) for stellar posi- 
tions, V magnitudes and spectral types (Spec). By 6 we denote 
the separations between the respective rosat cluster centre and 
star. 



BD 


HD 


#12000 


<5j2000 


e 


m v 


Spec. 


-00 301 


12134 


01:59:10.3 


+00:30:24 


1.94' 


8.28 


F0 


+ 18 315 


15551 


02:30:30.1 


+ 18:39:51 


3.59' 


8.25 


K0 


+28 1562 


67543 


08:09:34.3 


+28:11:46 


1.51' 


8.60 


F0 



solute value of the shear signal ruling out simultaneously higher 
values for C200 and r^. 

It will be interesting to see if the tentative trends in Fig. [9] 
will persist for a more complete cluster sample. We note that 
due to the substructure and non-sphericity of ACDM halos, we 
expect mass-concentration measured with weak lensing to be 
biased with respect to the more direct estimates from the sim- 
ulations (e.g., Bahe et al. 20121. At the current level of eight 



clusters, all mass estimates agree with each other within their 
error margins. We notice the relative statistical uncertainties to 
be only weakly dependent on the choice of C200 (free parameter, 
B01, or B12 mass-C2oo-relations). 

7.5. Masking of bright stars 

In Table [8] we summarise the properties of the magnitude 
8-9 stars that impede the analysis in the CL 0159+0030, 
CL0230+1836, and CL0809+2811 fields. By coincidence, 
these three most severe cases among the sample of 36 clusters 
are among our Megacam targets, reminding us that such fields 
must not be discarded when analysing a statistically complete 
sample. Using the example of CL 0809+28 11, we study the im- 
pact of these stars and their masking on the S -maps and mass 
estimates. 

Removing the masks generated for regions of deviant 
source density (Sect. |2.2| red squares in Fig. |D.3[ ) does not in- 
crease the number of usable galaxies significantly: Where scat- 
tered light strongly affects the local background estimation, 
sources are discarded in an early stage of catalogue preparation. 
At the position of the rosat centre, there are no detections in 
the first place. Without masking, the S -peak of CL 0809+28 11 
is shifted by 2'A to the north-east (closer to the rosat centre) and 
very slightly lower (5 max = 5.27 instead of S max = 5.39). Because 
our default model excises galaxies at < l'.5 separation, to avoid 
the strong lensing regime, the impact on the mass is below 2%. 

A more important point could be the extra uncertainty in 
the chosen cluster centre, as we find the largest offsets be- 
tween WL and X-ray peaks for clusters with large masks. The 
M 20o en < M 2oo Cen observed for CL 0159+0030 (Fig. [9} is likely 
due to a washed-out lensing peak. Nevertheless, we ao not ob- 
serve similar peculiarities for CL 0809+281 1. 

7.6. Further sources of uncertainty 

There are several potential sources of uncertainty which are not 



only for clusters imaged in g'r'i'. A reliable quantification of 
these errors will require further analysis. 

We also do not consider the uncertainty in the choice of 
max (|fi|) in the error analysis. However, we account for its effect 
via the shear calibration such that we do not expect a significant 
additional systematic error. Carefully calibrated simulations of 
cluster lensing are necessary to test our assumptions on the shear 
calibration factor. In Paper I, we observed in CL 0030+26 18 a 
counter-intuitive decrease of the best-fit value for r^j with in- 
creasing max(|e|). Indeed, only CL 1357+6232 shows a simi- 
lar relative decrease in r™ ( ". Averaging over all eight clusters, 
these cases are balanced by CL 0230+1836 and CL 1416+4446, 
for which we measure r™Q to increase with max(|e|). With 
max(|e|) = 1.0, we measure for four cases a smaller r^™ than 
for max (|e|) = 0.8, and in four cases a larger radius. The same 
holds for max(|e|) = 10 4 . These results suggest that the uncor- 
rected bias due to max (|e|) might be small. 

We notice that the roles of the shear calibration /o, con- 
sidered in Eq. ( 10 1 as cr ca ii and the correction f\{0) for cluster 
members cannot be completely disentangled. On the one hand, 
considering the cluster member correction separately is justified 
by the radial dependence of f\{9). On the other hand, we stress 
that the uncertainty in f\ (9) might be large due to the weak detec- 
tions of the cluster red sequence. In addition, the effect of cluster 
member decontamination on the mass estimate (3% to 7%; up 
to 11% using shear-peak centres) lies within the range of the 
related systematic error component cr~ ali which is significantly 
smaller than the statistical uncertainty in the mass. Hence, a pos- 
sible plan to consider cluster membership consistently - also for 
single-band clusters - would be to include it into the systematic 
error. Again, we suppose performing simulations of cluster WL 
fields to be helpful for the further investigations. 

Finally, we note that total lcr error intervals consistent with 
cluster masses close to zero do not mean these clusters are de- 
tected merely at the ~ lcr level: Shear calibration, ^-distribution, 
and triaxiality errors are multiplicative, such that they do not af- 
fect the detection significance. 



8. Summary and conclusion 

In this study, the second in the series on the 400d survey WL 
follow-up, we reduced and analysed MMT/Megacam observa- 
tions for seven clusters of galaxies. Building on Paper I, data 
reduction is performed using THELI, and WL shear catalogues 
are extracted using an implementation of the KSB+ algorithm. 

In the three cases, where we have MMT observations in 
g'r'i' , we define lensing catalogues based on a refined version 
of the three-colour method used in Paper I. By comparing with 
the colours observed for a Ilbert et aL| ( 2006[ > photo-z field, we 
exclude sources from regions in colour-colour-magnitude space 
containing a large fraction of foreground galaxies from the anal- 
ysis. For clusters with only one MMT band, we apply a cut in 
magnitude as background selection. 

We detect all of our 0.39 <z< 0.80 clusters using the aper- 
ture mass method (5 -statistics) at the >3.5<x level. Performing a 
Seitz & Schneider ( 200 1 ) l mass reconstruction, we find the pro- 



considered in Eq. ( 10 1, for instance, uncertainties in the determi- 



nation of the centres or the contamination correction available 



jected mass to follow the S -statistics closely. The WL masses of 
our clusters are determined from NFW modelling of their tan- 
gential shear profiles, yielding masses in the 10 14 Mq < M?L < 
2x 10 15 M© interval. 

Two of our clusters are exceptional due to their complicated 
shear morphology: For CL 1701+6414, where several known 
clusters lie close to one another in projection, we simultaneously 
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fit the shear of the two strongest S -peaks, identified with A 2246 
and our target. The field of CL 1641+4001 also exhibits multi- 
ple shear peaks, but we find no evidence for the presence of more 
than one cluster. 

An independent analysis of the CL 1701+6414 field us- 
ing archival CFHT data confirms the superposition of several 
weak lensing sources. By matching shear catalogues from MMT 
and CFHT, we find ellipticities measured with both instruments 
to be consistent with the assumption of noisy, but unbiased 
measurements of the same quantity. Hence, MMT/Megacam 
is proven to be equally good for WL science as the well- 
established CFHT/Megacam. We further produced a photo-z cat- 
alogue based on CFHT g'r'i'z' data of the field. Despite the shal- 
lowness of three bands, we are able to devise a coarse fore- 
ground/background selection for CL 1701+6414. This experi- 
ment again confirms not only shear peaks for several known 
clusters but also the validity of the magnitude-cut selection. 
Better photo-z data will be needed to potentially turn this cross- 
check into a calibration for single-band lensing data. 

We find the error budgets for our cluster masses to be dom- 
inated by statistical uncertainties (which can be suppressed by 
using a large cluster sample), but with a significant contribution 
of systematic uncertainties. Statistical uncertainties are naturally 
higher for relatively high-z clusters like ours, but the role of data 
quality (weather and instrumental conditions that led to dras- 
tic reductions in data depth) can hardly be overestimated. Our 
data show weak indications for clusters with low intrinsic con- 
centration to be more susceptible to model choices, such as the 
concentration parameter or assumed centre. 

As the second paper in the 400d WL series, the main results 
of this study are that: 

- Instrumental effects are well under control. 

- Reliable masses can be obtained in the presence of bright 
stars close to the cluster centre - an important finding for a 
successful follow-up of a complete sample! 

- Ground-based WL works at least till z»0.8. 

- We can correct for extreme cases of massive foreground 
structures. 

- Clusters of low concentration are a common occurrence in 
the high redshift, moderate mass population. 

Furthermore, we identify areas of possible future improvements 
of the methods we applied: As observational constraints will for- 
bid complete homogeneity of the data analysis to some degree 
for each survey. Therefore, better methods to calibrate WL anal- 
yses in particular with a different number of available filters need 
to be developed. This applies specifically to background selec- 
tion and correction for cluster members. 

Concerning the modelling of clusters with complex shear 
morphology, we point out that low-mass clusters with a strong 
fraction of disturbed or merging systems will build the bulk of 
the population observed by deep and wide future surveys like 
eRosita or Euclid. Hence, the question arises how clusters de- 
viating from a simple NFW mass distribution can be weighed 
most accurately. A possible method is to apply methods that 
do not assume radial symmetry, e.g. aperture mass techniques. 
Nevertheless, profile-fitting methods are well established and in- 
creasingly well understood. Alternatively, by combining simu- 
lational efforts with improved data analysis and modelling, bi- 
ases resulting from profile assumptions can be corrected. Thus, 
applying profile fits even to low-signal, merging clusters might 
prove the best method to measure reliable masses for large clus- 
ter samples. 



Acknowledgements. HI likes to thank Frank Bertoldi for support of this 
work and Matthias Klein, Reiko Nakajima, Mischa Schirmer, Ismael Tereno, 
Bharadwaj Vijaysarathy, Daniela Wuttke, and Yu-Ying Zhang for helpful 
discussions. Partial support for this work has come from the Deutsche 
Forschungsgemeinschaft (DFG) through Transregional Collaborative Research 
Centre TRR 33 as well as through the Schwerpunkt Program 1 177. THR ac- 
knowledges support from the DFG through Ffeisenberg grant RE 1462/5 and 
grant RE 1462/6. CLS was supported in part by Chandra grants GO9-0135X, 
GO9-0148X, and G01-12169X, and by the F. H. Levinson Fund of the Silicon 
Valley Community Foundation, which helped support the MMT observations. 
We acknowledge the grant of MMT observation time (program 2007B-0046) 
through NOAO public access. 



References 

Adelman-McCarthy, J. K., Agiieros, M. A., Allam, S. S., et al. 2008, ApJS, 175, 
297 

Allen, S. W., Evrard, A. E., & Mantz, A. B. 201 1, ArXiv astro-ph/1 103.4829 

Bahe, Y. M., McCarthy, I. G, & King, L. J. 2012, MNRAS, 421, 1073 

Barkhouse, W. A., Green, P. J., Vikhlinin, A., et al. 2006, ApJ, 645, 955 

Bartelmann, M. 1996, A&A, 313, 697 

Benftez, N. 2000, ApJ, 536, 571 

Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 

Bhattacharya, S., Habib, S., & Heitmann, K. 2011, ArXiv astro- 

ph.CO/1 112.5479 (B12) 
Blanchard, A. 2010, A&A Rev., 18, 595 
Bohringer, H. & Werner, N. 2010, A&A Rev., 18, 127 
Bolzonella, M., Miralles, J., & Pello, R. 2000, A&A, 363, 476 
Bullock, J. S., Kolatt, T. S., Sigad, Y, et al. 2001, MNRAS, 321, 559 (B01) 
Burenin, R. A., Vikhlinin, A., Hornstrup, A., et al. 2007, ApJS, 172, 561 
Burenin, R. A. & Vikhlinin, A. A. 2012, ArXiv astro-ph.CO/1202.2889 
Cohn, J. D. & White, M. 2005, Astroparticle Physics, 24, 316 
Coleman, G. D., Wu, C.-C, & Weedman, D. W. 1980, ApJS, 43, 393 
Corless, V. L. & King, L. J. 2007, MNRAS, 380, 149 
Corless, V. L. & King, L. J. 2009, MNRAS, 396, 315 

Dietrich, J. P., Bohnert, A., Lombardi, M., Hilbert, S., & Hartlap, J. 201 1, ArXiv 

astro-ph/1 103.4607 
Dietrich, J. P., Erben, T., Lamer, G, et al. 2007, A&A, 470, 821 
Donahue, M., Scharf, C. A., Mack, J., et al. 2002, ApJ, 569, 689 
Eke, V. R., Cole, S., & Frenk, C. S. 1996, MNRAS, 282, 263 
Erben, T, Hildebrandt, H., Lerchster, M., et al. 2009, A&A, 493, 1197 
Erben, T, Schirmer, M., Dietrich, J. P., et al. 2005, Astronomische Nachrichten, 

326, 432 

Erben, T, van Waerbeke, L., Bertin, E., Mellier, Y, & Schneider, P. 2001, A&A, 
366,717 

Goto, T, Sekiguchi, M., Nichol, R. C, et al. 2002, AJ, 123, 1807 
Hartlap, J., Schrabback, T., Simon, P., & Schneider, P. 2009, A&A, 504, 689 
Hetterscheidt, M., Erben, T., Schneider, P., et al. 2005, A&A, 442, 43 
Heymans, C, van Waerbeke, L., Bacon, D., et al. 2006, MNRAS, 368, 1323 
High, F. W., Hoekstra, H, Leethochawalit, N., et al. 2012, ArXiv e-prints 
High, F. W., Stubbs, C. W., Rest, A., Stalder, B., & Challis, P. 2009, AJ, 138, 110 
Hildebrandt, H., Erben, T., Dietrich, J. P., et al. 2006, A&A, 452, 1121 
Hoekstra, H. 2003, MNRAS, 339, 1155 

Ilbert, O., Arnouts, S., McCracken, H. J., et al. 2006, A&A, 457, 841 
Ilbert, O., Capak, P., Salvato, M., et al. 2009, ApJ, 690, 1236 
Israel, H., Erben, T, Reiprich, T. H., et al. 2010, A&A, 520, A58 
Kaiser, N. & Squires, G. 1993, ApJ, 404, 441 
Kaiser, N., Squires, G, & Broadhurst, T. 1995, ApJ, 449, 460 
Kasun, S. F. & Evrard, A. E. 2005, ApJ, 629, 781 
Komatsu, E., Smith, K. M., Dunkley, J., et al. 201 1, ApJS, 192, 18 
Kowalski, M., Rubin, D., Aldering, G, et al. 2008, ApJ, 686, 749 
Kravtsov, A. V, Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 
Larson, D., Dunkley, J., Hinshaw, G, et al. 2011, ApJS, 192, 16 
Lopes, P. A. A., de Carvalho, R. R., Gal, R. R., et al. 2004, AJ, 128, 1017 
Mantz, A., Allen, S. W, Rapetti, D., & Ebeling, H. 2010, MNRAS, 406, 1759 
McLeod, B. A., Conroy, M., Gauron, T. M., Geary, J. C, & Ordway, M. P. 2000, 
in Further Developments in Scientific Optical Imaging, ed. M. B. Denton, 1 1 
Medezinski, E., Broadhurst, T., Umetsu, K, et al. 2010, MNRAS, 405, 257 
Meneghetti, M., Rasia, E., Merten, J., et al. 2010, A&A, 514, A93 
Miller, C. J., Nichol, R. C, Reichart, D., et al. 2005, AJ, 130, 968 
Nagai, D., Kravtsov, A. V, & Vikhlinin, A. 2007a, ApJ, 668, 1 
Nagai, D., Vikhlinin, A., & Kravtsov, A. V. 2007b, ApJ, 655, 98 
Navarro, J. F, Frenk, C. S., & White, S. D. M. 1995, MNRAS, 275, 720 
Navarro, J. F, Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 
Navarro, J. F, Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Oguri, M., Takada, M., Okabe, N., & Smith, G. P. 2010, MNRAS, 405, 2215 
Pillepich, A., Porciani, C, & Hahn, O. 2010, MNRAS, 402, 191 



21 



Holger Israel et al.: The 400d weak lensing survey II 



Pillepich, A., Porciani, C, & Reiprich, T. H. 201 1, ArXiv astro-ph/1 1 1 1.6587 
Planck Collaboration, Aghanim, N, Arnaud, M., et al. 2012, ArXiv astro- 
ph.CO: 1204.2743 

Plionis, M., Basilakos, S., Georgantopoulos, I., & Georgakakis, A. 2005, ApJ, 
622, L17 

Predehl, P., Andritschke, R., Bohringer, H., et al. 2010, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7732, 
Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series 

Press, W. H. & Schechter, P. 1974, ApJ, 187, 425 

Rasia, E., Ettori, S., Moscardini, L., et al. 2006, MNRAS, 369, 2013 

Reimers, D., Toussaint, E, Hagen, H., Hippelein, H., & Meisenheimer, K. 1997, 

A&A, 326, 489 
Reiprich, T. H. 2006, A&A, 453, L39 

Romer, A. K., Nichol, R. C, Holden, B. P., et al. 2000, ApJS, 126, 209 
Rosati, P., Borgani, S., & Norman, C. 2002, ARA&A, 40, 539 
Sarazin, C. L. 1988, X-ray emission from clusters of galaxies (Cambridge 
University Press) 

Schirmer, M., Erben, T., Hetterscheidt, M., & Schneider, P. 2007, A&A, 462, 
875 

Schneider, P. 1996, MNRAS, 283, 837 

Schneider, P. 2006, "Weak Gravitational Lensing" in: Gravitational Lensing: 
Strong, Weak and Micro: Saas-Fee Advanced Courses, Volume 33 (Springer- 
Verlag Berlin Heidelberg), 269 ff. 
Schneider, P., King, L., & Erben, T. 2000, A&A, 353, 41 
Schrabback, T., Erben, T., Simon, P., et al. 2007, A&A, 468, 823 
Schrabback, T., Hartlap, J., Joachimi, B., et al. 2010, A&A, 516, A63 
Schuecker, P. 2005, in Reviews in Modern Astronomy, Vol. 18, Reviews in 

Modern Astronomy, ed. S. Roser, 76-105 
Seitz, S. & Schneider, P. 1996, A&A, 305, 383 
Seitz, S. & Schneider, P. 2001, A&A, 374, 740 

Stanek, R., Rasia, E., Evrard, A. E., Pearce, E, & Gazzola, L. 2010, ApJ, 715, 
1508 

van Waerbeke, L., Mellier, Y., Radovich, M., et al. 2001, A&A, 374, 757 
Vikhlinin, A., Burenin, R. A., Ebeling, H., et al. 2009a, ApJ, 692, 1033 
Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009b, ApJ, 692, 1060 
Vikhlinin, A., McNamara, B. R„ Forman, W., et al. 1998, ApJ, 502, 558 
Voit, G. M. 2005, Reviews of Modern Physics, 77, 207 

von der Linden, A., Best, P. N., Kauffmann, G., & White, S. D. M. 2007, 

MNRAS, 379, 867 
Wright, C. O. & Brainerd, T. G. 2000, ApJ, 534, 34 
Zhang, Y.-Y, Finoguenov, A., Bohringer, H., et al. 2008, A&A, 482, 451 
Zhang, Y.-Y, Okabe, N, Finoguenov, A., et al. 2010, ApJ, 711, 1033 



Appendix A: Photometric calibration details 



the fit parameters have been determined, applying an iterative 



3<r-clipping fit of Eq. (A.l 



Applying the Hildebrandt et al. ( 2006 ) method, photometric cal- 
ibration of our data is established by fitting instrumental (mi nst ) 
to reference (wsdss) magnitudes for a sample of objects, taking 
into account variable airmass a and a colour term csdss describ- 
ing the transformation between Megacam and SDSS filter sys- 
tems: 

Winst - '«SDSS = y8fCSDSS,f + JiO + Zf , (A. 1) 

Depending on the photometric quality of the observations, we fit 
the zeropoint Zf together with the parameters fit or jf in optimal 
conditions, or keep jf fixed at the default value depending on the 
filter / for poorer conditions. 

The resulting values for the fit parameters, as well as the 
colou r ind ices csDSS.f for the different filters are presented in 
Table 



A.l 



for all photometric nights of our Megacam runs p] We 
find the zeropoints Z t of the photometric nights to agree among 
the g'r'i' filters, with a largest deviation of «0.15 mag. The scat- 
ter Am' = mjnst - ;«sdss + ySfCsDSS.f + JfCi + Zf of the individual 

i has a com- 



SDSS standards about the best-fit solution (Fig. A.l 
parable amplitude. The errors of Zf given in Table 
formal fitting errors. Figure IaT| presents the data from which 



A.l are the 



11 Note that some values in Table |AjJ are corrected w.r.t. Table A.l in 
Paper I. The amount of these corrections is of the order of, and in most 



Comparing the colour terms /Jf for the different nights, we 
find considerable agreement within the values for each of the 
three bands, although the formal errors underestimate the true 
uncertainties. We suggest that the large span in values of j3 g 
might be caused by the known dependence of the filter through- 
put on the distance to the optical axis. Plotting the scatter Am' 
as a function of the separation 6q of the source from the optical 
axis of Megacam (Fig. |A. 1 1 ), given by the pointing position in the 
fits header, we can confirm trends of Am'(#o) in all filters, most 
pronounced for the g' band data taken on 2005 October 30. This 
trend is likely caused by a combination of the sky concentra- 
tion effect (position-dependent illumination due to scattering in 
the telescope optics) and the position-dependent transmissivity 
of the Megacam filters, which is strongest in the g' band (cf. Fig. 
A. 3 in Paper I). A more conclusive investigation of this issue, 
requiring full propagation of errors on instrumental magnitude, 
lies beyond the scope of this paper. Because the radial depen- 
dence observed in Fig. 



A. 1 does not exceed the residual scatter 



for sources at the same 6q, the global photometric fits (Eq. A.l 
fulfil the requirements of our analysis. 



Appendix B: Details of Background Selection 

In the intermediate magnitude range mbright <r'< mf a ; nt , our lens- 
ing catalogues for three-band clusters include galaxies selected 
from g'-r' versus r'-i' colour-colour- diagr ams (Sect.|2.4|l. We 
find our method justified by considering the |Ilbert et al.| ([2006 ) 
photo-z catalogue: 

Figure [B~T| presents the galaxy numbers and the fraction of 
z p h < 0.50 sources in the Deep 1 photo-z catalogue as a func- 
tion of the r' magnitude and g'-r' and r'-i' colours. First, the 
catalogue is divided into its dodeciles in r', i.e. twelve magni- 
tude bins of equal population are defined where the k-th bin con- 
sists of the galaxies r' k x < r' < r' k . By r' k , we denote the mag- 
nitude of a source such there is a fraction of k/12 of brighter 
galaxies in the catalogue. Second, for each dodecile, we show 
the number Ntj of galaxies falling into grid cells of mesh size 
A(g' -r') = A(r' — £') =0.1, using a grey scale. Figure B.l high- 
lights that at bright r', only a narrow strip in the colour-colour 
space spanned by g' - r' and r' - i' is populated, while the lo- 
cus of galaxies becomes much more diffuse towards fainter r'. 
Third, for each grid cell, we determine the fraction of galax- 
ies we define as foreground sources, i.e. the sources with a red- 
shift estimate z p h < 0.40. The red, green, and blue contours in 
Fig. B.l mark regions of the colour-colour space populated by 
25%, 50%, and 75% of foreground galaxies compared to the 
z d = 0.40 clusters, CL 0159+0030 and CL 0809+2811. The con- 



tours are defined such that /f g exceeds the respective threshold 
in all grid cells enclosed by the contour. 

As expected, /f g generally decreases towards fainter mag- 
nitudes, with only a few z p h < 0.40 sources at r' > 26.0. For 
all magnitudes, foreground sources with r' - i' > 0.5 are rare. In 
the brightest three dodeciles, a well-defined region with a dis- 
tinctive edge towards redder r'-i' colours exists^] Although the 
preferred locus of z p h < 0.40 galaxies depends little on the r' 
magnitude, the zone populated by low-z objects becomes more 



12 Towards very blue r' - i' colours, few galaxies are found in the 
CFHTLS D 1 catalogue, basically all of them at low z < 0.4 redshift. This 



cases smaller than, the scatter observed in Fig. A.l 



can be seen from the contours in Fig. |B.1| which follow the irregular 
shape of the point cloud. We choose a conservative min(r'-;') = -1.0 
limit for the selection polygons. 
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Fig. A.l. Accuracy of the photometric calibration: For the different combinations of filters and nights used to calibrate the data sets 
discussed in this work and Paper I, the scatter Am' around the best-fit solution (solid line) is shown. Each point corresponds to an 
SDSS standard source for which the abscissae give the separation 9$ in arc minutes from the centre of the pointing. Note that for 
each panel a maximum Am' has been determined by iterative 3<x-clipping. 
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Table A.l. Coefficients of photometric calibration defined by Eq. (A.li for all photometric nights within our 
MMT/Megacam 400d observations. 



Filter Obs. Date 






A 


CSDSS 


7f 




g' 2005-10-30 


27.277 ± 


0.005 


0.106 ± 


0.007 


g'-r' 


(-0.15) 8 


2 


2005-11-01 


27.286 ± 


0.005 


0.116 ± 


0.005 


g'-r' 


(-0.15) § 


2 


V 2005-10-31 


26.426 ± 


0.002 


0.124 + 


0.002 


r' - i' 


(-0.05)8 


2 


2005-11-01 


27.408 ± 


0.009 


0.119 ± 


0.002 


r' - i' 


-0.03 ± 0.01 


3 


r 1 2005-06-07 


26.819 ± 


0.001 


0.040 ± 


0.001 


g'-i' 


(-0.10)8 


2 


2005-06-08 


26.834 ± 


0.008 


0.048 ± 


0.001 


g'-i' 


-0.12 ±0.01 


3 


2005-10-30 


26.950 ± 


0.018 


0.046 ± 


0.002 


g'-i' 


-0.10 ±0.02 


3 


2005-10-31 


26.959 ± 


0.004 


0.042 ± 


0.003 


g'-i' 


(-0.10)8 


2 


2005-11-01 


26.960 ± 


0.008 


0.048 ± 


0.004 


g'-i' 


(-0.10)8 


2 


2005-11-08 


26.807 ± 


0.005 


0.046 ± 


0.003 


g'-i' 


(-0.10)8 


2 



1 Normalised to an exposure time of Is and an airmass a = 0. 
* Number of parameters used in the fit. 
8 Fixed at the default value. 

Table B.l. Cuts defining the polygons used for background selection for the z ~ 0.40 and z = 0.80 clusters, based on the colours of 
foreground galaxies (Fig. |B.l| . We specify the values of g'-r, r'-i', and sp=/3(r'-i') - (g'-r') at the edges of the exclusion polygons 

formb r ight<r'<m fa i n t. 



Redshift 


min(r'-i') max (r'-i') 


min (g'-r') 


max (g'-r') 


yS 


min(.s /3 ) 


max (sp) 


z = 0.4 


-1.0 0.7 


0.1 


1.6 


2.5 


-3.5 


0.5 


z = 0.8 


-1.0 1.2 


0.3 


1.7 


1.5 


-3.0 


0.7 



B.l 



noticing 



diffuse for fainter sources. The insignificant role of foreground 
galaxies r' > 25.0 justifies that our background selection includes 
all galaxies fainter than mf a ; nt . Furthermore, the secondary role of 
wibright compared to mf a ; nt becomes clear from Fig. 
the small number of r' <20 galaxies. 

Calculating ff s for a cluster redshift of Zd =0.80, the regions 
in colour-colour space where a given value of /f g is exceeded 
extend as well towards fainter r' as towards redder g' - r' and 
r' - i' colours. For z& - 0.80, only a small number of background 
sources remain in the first dodecile, while there are significant 
foreground objects even in the r' w < r' < r' x , bin. 

We adjust our background selection polygon to Zd = 0.40 
and za = 0.80 by defining criteria based on /f g > 0.25 contours in 
the three brightest dodeciles (r' < 23.91) of the Deep 1 photo-z 
catalogue (Table B.l and yellow polygon in Fig. B.l i. We ex- 



clude galaxies at m^ght ^ f' < mf a ; nt and matching these crite- 
ria from the lensing catalogues. Performing a cross-check for 
Zd = 0.50, we confirm the background selection in Paper I to be 
sensible, although not optimal. In fact, more "self-calibrations" 
can be achieved by combining three-colour photometry with 
photo-z catalogues (Klein et al. in prep.). 



Appendix C: Details of photo-z analysis 

C.1. Spectroscopic calibration 

Because our photo-z catalogue for the CFHT CL 1701+6414 
field was distilled from only four bands, three of which have 
rather shallow exposure time (Table [5}, we tested its quality by 
comparison with publicly available SDSS spectroscopy redshifts 
of the same field. SDSS spectra are only available for a selection 
of the brightest (r' < 19) sources, with a total of 88 matches for 
the as 270000 object photo-z catalogue. Out of the 58 sources 
flagged as good by BPZ in all four filters, 45 are identified as 
normal galaxies by SDSS, nine are identified as stars, and four 
as QSOs. Figure C.l displays z p h as a function of z spec of the 



matched sources for which the BPZ quality parameter ODDS is 
o > 0.7. The size of the symbols in Fig. C.l (star symbols for 
stars, filled circles for normal galaxies, and filled squares for 
QSOs) corresponds to the value of o. 

Generally, the photo-z uncertainties (given by the BPZ 
Z_B_MIN and Z_B_MAX parameters) are large, although the 



sources in Fig. C.l rank among the brightest in the catalogue. 
Nevertheless, the z p h estimates for normal galaxies seem to fol- 
low a remarkably narrow and mono tonic function of z speC i in the 
range < z spe c 55 0.45 probed by the SDSS spectral targets in 
Fig. C.l The step at z spe c ~ 0.10, below which galaxies get as- 
signecTzph ~ 0.0 and above which they are overestimated to be at 
z p h ~ 0.3, can be explained by the lack of a u' filter crucial for 
detecting the 400 nm break at these redshifts. In particular, this 
applies to galaxies in the z « 0.22 structures in the foreground to 
CL 1701+6414, as exemplified by the z p h ~ 0.40 for one of the 
A 2246 BCG candidates at z spe c -0.235. 

Photo-z estimates for the highest redshift (z ~ 0.45) galaxies 
with SDSS spectra are stunningly accurate, despite the large un- 
certainties. The overall trend seen in Fig. C.l is consistent with 
the results of the CFHTLS-Archive-Research Survey ( |Erben| 
et al. 2009). In their analoguous comparison of CFHT BPZ 
photo-zs to SDSS spectra, they find a turnover to z p h < z spec for 
Zspec ^0.45 (for a small absolute number of such galaxies). 

It is a lucky coincidence that the z p h-z spe c-relation inter- 
sects the dotted equality line precisely at the redshift of our 
cluster of interest. There are six galaxies 0.448 < z spe c < 0.457 
among the SDSS spectral targets (inlay in Fig. C.l i, while there 
are none in the 0.36 < z spe c < 0.44 range. These six include the 
BCG of CL 1701+6414 at z spec = 0.4523 ± 0.0001, for which 
BPZ returns z p h = 0.45 ± 0.19. However, even the closest of the 
other five is separated by 12'0 or 4.0 Mpc in projection and thus 
not part of or closely interacting with CL 1701+6414. Still, its 
Zspec = 0.4522 ± 0.0001 indicates they might belong to the same 
large-scale structure. 
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Fig. B.l. Fraction of z p h < 0.40 galaxies in the Deep 1 (Ilbert et al. | |2006| l field as a function of their g' — r' and r' — i' colours 
and r' magnitude. Each panel shows a dodecile of the photo-z catalogue, i.e. one of twelve equally populated magnitude bins, 
where the £-th dodecile includes all galaxies r' k l < f < r' k . In each panel, the number iVy of galaxies within cells of mesh size 
A(g' — r') = A(r' — V) = 0.1 is shown, using the same grey scale oc ^JN~ij. White grid cells are empty. The red, green, and blue contours 
enclose regions in which 25% (50%, 75%) of galaxies have a z p h < 0.40. Based on the distribution of z p h<0.40 galaxies in the three 
brightest dodeciles, we define the yellow polygon (see Table B.l i, in order to remove foreground galaxies from the CL 0159+0030 
and CL 0809+28 11 fields. 



The photo-z estimates for QSOs (right panel of Fig. C.l i are 
clearly off, which doesn't come as a surprise as their is no QSO 
spectrum among the templates employed by BPZ. Although five 



out of eight stars in Fig. C. 1 (left panel) get assigned z p h < 0.02, 
there are also three cases in which our photo-z catalogue sug- 
gests higher z p h for objects that by visual inspection and SDSS 
classification clearly stellar: The available CFHT photometry 
alone does not allow for an accurate star/galaxy classification. 



C.2. Matching with MMT 

To investigate what benefit the four-band photo-zs yield, once 
the usual selection of galaxies by magnitude and half-light radius 



(Sect. |2.4| i is applied, we now match the photo-z catalogue with 
the MMT KSB catalogue. We first notice that although 95.6% 
of galaxies in the MMT KSB catalogue are matched bijectively 
to a CFHT photo-z source, only 20.4% of the MMT lensing 
sources satisfy the quality criteria of four usable CFHT bands 
(NBPZ_G00DFILT= 4), o > 0.8, and no masking (CANDMASK= 0). 
We call these estimates high-quality photo-zs. 

Counter-intuitive at the first glance, 76% of the high-quality 
matches have z p h > 0.5, i.e. are likely background galaxies to 
CL 1701+6414. However, this effect can be traced back to the 
near-absence of high-quality estimates of z p h < 0.3. This is con- 
sistent with our expectation from the comparison to the spectro- 
scopic redshifts (Fig. C.l I, where the lack of a u' filter systemat- 
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Fig.C.l. CFHT photometric red- 
shifts plotted against spectroscopic 
redshifts from SDSS. The left 
panel (star symbols) shows the 
z p h for objects identified as stars 
in SDSS, the middle panel (filled 
j= circles) for normal galaxies, and 
N the right panel (filled squares) 

for QSOs. A small inlay gives a 
zoomed version for z spi;c ~ 0.45 
galaxies. The size of the symbols 
marking the photo-z estimate cor- 
respond to the quality parameter o 
(ODDS): big symbols for o > 0.9, 
medium-sized symbols for 0.8 < 
o < 0.9, and small symbols for 
0.7 < o < 0.8. Error bars for 
the Zph give the uncertainty inter- 
val defined by the BPZ Z_B_MIN 
and Z_B_MAX parameters. Objects 
inside a flagged region of one of 
the CFHT images (CANDMASK= 1) 
are shown in grey. 
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Fig. C.2. Upper panel: Histograms of the MMT-photo-z 
matched catalogue as a function of MMT r' magnitude: Plotted 
are all matches (thick black line) and the subset of matches with 
photo-zs passing all quality criteria ("high quality", blue line). 
Lower panel: Fractions of galaxies z p h > 0.45 within: all MMT- 
photo-z matches (thick black line), the "high quality"matches 
(blue line), and of high-quality matches z p h > 0.45 within all 
(dashed red line). A thin dotted line gives the magnitude cut at 

'«faint = 21.9. 



containing brighter galaxies on average, traces a different pop- 
ulation than our MMT lensing catalogue. 



The upper panel of Fig. C.2 showing the magnitude distri- 
butions of the high-quality photo-z catalogue (solid blue line) 
compared to all matches (thick black line) demonstrates that 
high-quality photo-zs tend to belong to brighter galaxies. This 
can be seen from the modes of the histograms and is not surpris- 
ing given the necessary detection in the shallow g'i'z' images. 
While 98.5% of the sources in the MMT galaxy catalogue get 
matched to a photo-z galaxy, only for 26.7% the photo-z passes 
all quality cuts. 

In particular, the decline with magnitude of the fraction of 
high-quality matches affects the fraction of background galax- 
ies with respect to CL 1701+6414 at z = 0.45 using our photo- 
zs: Considering high-quality matches only (solid blue curve in 
the lower panel of Fig. C.2 1, the fraction of "photometric back- 
ground" (z p h >0.45) increases strongly with MMT r' magnitude. 
In fact, all of the few r' MMr > 25.8 high-quality matches show 
Zph > 0.45. With respect to the complete catalogue (thick black 
line), however, the fraction of photometric background galaxies 
peaks at r^ MT w 24 and ~ 0.6 and decreases towards fainter r' MMJ . 
The fraction of high-quality z p h > 0.45 galaxies compared to all 
matches (red dashed line in Fig. C.2 1 runs rather flat with r' MMJ , 



never exceeding 0.3 and subsuming only 15.9% of all matches. 

We conclude that the quality of the CFHT data entering the 
photo-z estimation makes possible a rough estimation of a nor- 
mal galaxy's redshift, i.e. to decide if it is more likely to be in the 
foreground or in the background, but not a precise redshift dis- 
tribiution from which (fi) could be inferred more precisely than 
using a proxy photo-z catalogue of high quality (Sect.|23]l. 



ically offsets z p h for z spe c < 0.45 galaxies to higher values. The z p h 
distribution derived from the four-band CFHT data deviates far 
from the redshift distributions known from well-studied photo-z 
fields ( jllbert et al.|2006l[2009] also see Fig. B.6 in Paper I). This 
holds in particular for the high-quality sub-catalogue which, 



C.3. A photo-z shear catalogue 

The photo-zs drawn from the CFHT g'r'i'z' bands provide 
us with a rough redshift estimate. In order to test whether 
this information can be used to disentangle the shear signals 
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Fig.C.3. 5-statisics in the CL 1701+6414 
field using the CFHT photo-z catalogue. 
Thick black contours mark the shear sig- 
nal from the z p h < 0.45 catalogue, medium- 
thick blue contours are derived from the 
complementary z p h ^ 0.45 sources. Thin 
orange contours show the signal from the 
complete MMT galaxy catalogue. All con- 
tours start at S = 1 and are spaced by 



AS = 1. The underlying image and other 
contours are the same as in Fig. [3] 



of CL 1701+6414 and the foreground structure, in particular 
A 2246, we divide the MMT galaxy catalogue: Galaxies with 
z p h < 0.45 are sorted into the "photo-z foreground" catalogue, 
galaxies with z p h > 0.45 are sorted into the "photo-z background" 
catalogue. Because of the poor quality of most photo-z estimates 
we expect only a crude selection. 

Figure [C3] shows the S -maps resulting from these two cat- 
alogues, overlaid on the MMT r'-image in the same fashion 
as for Fig. [3] Thick black contours denote iso-S -contours from 
the photo-z background catalogue, including 49.9% of the lens- 
ing catalogue (10.3 galaxies/arcmin" 2 ). Solid blue contours in 
Fig. |C.3| are drawn from the complementary photometric fore- 
ground catalogue; the signal from the complete galaxy catalogue 
is shown as thin orange contours. 

The morphology of the S -peaks in the photo-z background 
map follows in its main features the complete catalogue, as we 
expect from a sample of true z>0.45 galaxies. With S max = 3.14, 
the peak to be associated with CL 1701+6414 is nearly as strong 
as for the full catalogue, and closer to the cluster's rosat posi- 
tion. The A 2246 peak shows a similar high fraction of the com- 
plete catalogue signal, but the two clusters appear to be better 
separated. The photo-z-foreground S -morphology bears little re- 
semblance to Fig. [3] Although we still measure S ~ 2 close to the 
position of CL 1701+6414, it can not be seen as a distinct peak. 
A 2246 is detected just below 3cr, with S > 3 only measured 
for the "SW peak". This is consistent with our expectations: As 
A 2246 is at lower redshift, some signal should persist in a true 
foreground catalogue. 

Re-defining the photo-z catalogue such that it only con- 
tains high-quality photo-zs > 0.45 (cf. Sect. |C.2[) results in a 



good resolution between the S -signals of the two main clus- 
ters but such catalogue suffers from the sparsity of sources (4.8 
galaxies/arcmkr 2 ). 

Despite the outcome of this experiment matching our ex- 
pectations, we keep in mind the typical uncertainty of cr(z ph ) » 
0.25 even for the high-quality photo-zs (Fig. C.l i, similar to the 



redshift separation of CL 1701+6414 and A 2246. Hence, the 
photo-z selection using the available data is not inherently bet- 
than the Sect. 



ter 



4.1 



magnitude cut. Nonetheless, the CFHT 
photo-zs and lensing measurements confirm the detection of 
CL 1701+6414 as a shear source distinct from A 2246 and give 
credibility to its mass estimate, the aim of our investigations. 



Appendix D: Notes on individual clusters 

D.1. CL 01 59+0030 

Being located in the SPSS equator ial strip, CL 0159+0030 has 
been detected by |Goto et al.| ( |2002] > in the SDSS commissioning 
data, using their photometric "cut and enhance" cluster finder. 
Plionis et al. (2005) followed up |Goto et al. ( 2002| l cluster can- 
didates using archival XMM-Newton observations. From the 
3800 s PN observation Plionis et al. (2005) analysed, only a 3cr 



upper flux limit of 2. 1 X 10" 14 erg cm~ z s _1 in the 0.5-2.0 keV en- 
ergy range could be inferred. This non-detection disagrees both 
with th e flux of 3.3 ± 0.4 x 10° erg cm" 2 s "' Vikh linin et al.| 
(2009a) measure for CL 0159+0030 with rosat and with their 
chandra flux of 3.6 x 10 13 erg cm" 2 s" 1 in the same band. 
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D.2. CL 0230+1 836 

For CL 0230+ 1836, there are neither detections of the cluster it- 
self, independent from the 400d survey, nor other galaxy clusters 
within a 20' radius listed in NED. To our knowledge, we are the 
first to study this high-z cluster with deep optical observations. 

D.3. CL 0809+2811 

We hypothesise that CL 0809+28 11 is identical to 
ZwCl 0806.5+2822 at a J20 oo = 08 h 09 m 34 s , 6 mm = +28° 13a, 
a position 1'.9 off the CL 0809+28 11 rosat centre and at 
similar distance to the bright star in the field, where we do 
not see a concentration of galaxies. Neither do we observe an 
overdensity of galaxie s at the position of a secondary shear 
peak with S = 2.9 (Fig. |D3]i. It is located at a J2 ooo = 08 h 09 m 08 s , 
<5j2000 = +28 °05'22". No cluster within 3' of this position is 
known to NED. 



D.4. CL 1357+6232 



Lopes et al. ( 2004) > conducted a cluster survey on digitised 
Second Palomar Observatory Sky Survey plates, using a Voronoi 
tesselation technique. In their catalogue, they quote a clus- 
ter of galaxies at a J20 oo = 13 h 57 m 22 s , 6 s2 ooo = +62°33'11", 
where there is no source in the Megacam image. Using the 
relation found between r' magnitude, g' - r' colour an d z spec 
for a subsample of clusters with spectroscopic redshifts, [Lopes 



et al.| ( |2004| l assign z = 0.19 to their detection. (NSCS 
J 135722+6233 11, their #7243). Noting that the position of 
NSCS J 135722+6233 11 is only 16" from the rosat centre of 
CL 1357+6232, we speculate that it might be the result of a con- 
fusion of CL 1357+6232 with two bright galaxies to its east, one 
of which (SDSS J135723.83+623246.1) has a measured redshift 
of z = 0.078. 



the strongest lensing detection, but also appears to be the op- 
tically richest system in the Megacam r'-band image. Therefore, 
CL 1416+4446 possibly presents the most massive system in a 
physically interacting super-structure, indicated by the z ~ 0.40 
redshifts of all mentioned clusters. Judging by the X-ray mor- 



phology, Vikhlinin et al. ( 2009a I classify CL 1416+4446 as a re 



laxed system, though. We conclude that the CL 1416+4446 field 
qualifies as an interesting candidate for further investigation con- 
cerning the existence of a super-cluster. 

D.6. Cross-identifications of CL 1 701 +641 4 



Vikh linin et al.| (1998i, on whose 160d rosat catalogue the 
400d sample builds (cf. Burenin et al.|2 007 ) detect four clusters 
in the field: VMF 189 (A 2246), VMF 190 (CL 1701+6414), 
VMF 191, and VMF 192. In an independent rosat analysis, 
Donahue et al.| ( |20Q2| > detect these same four clusters plus 
RX J 1702+6407, which we do not detect in WL. The redshift of 



z = 0.7 found for the |Donahue et al. (2002) optical counterpart of 
CL 1701+6414 deviates from the redshift of z — 0.45 measured 



by |Burenin et al.| ( |2007| l and all other references. CL 1701+6414 
is further listed as RX J1701 .3+6414 in the Bright Serendipitous 
High-Redshift Archival rosat Cluster sample (Bright SHARC, 
|Romer et al.|2000) . 



D.5. CL 1416+4446 

In addition to CL 1416+4446, we detect two other shear peaks at 
> 3cr significance to the west and south-west of CL 1416+4446. 
|Lopes etal.|p004) list a cluster NSCS J141623+444558 in their 



catalogue wh ich, by NED , is id entified with CL 1416+4446. 
Furthermore, Lope s et al. (2004) detected a cluster of galax- 



ies at a J20 oo = 14 h 16 m 09 s and feooo = +44°38'51", with a 
redshift of z — 0.39. Less than 2' north-east of these coor- 
dinates we find the south-western shear peak which coincides 
with the g' = 2 0.1 galaxy SDSS J141613.33+44 395 1.3. For this 
source, SDSS (Adelman-McCarthy et al. 2008) quotes a spec- 
troscopic redshift of z — 0.397. Note that the brighter galaxy 
SDSSJ141603.01+443725.1, located 2' further to the south- 



west from the Lopes et al. ( 2004)1 cluster position has an SDSS 



Zspec =0.310 and does not correspond to an S -peak. 



Barkhouse et al. (2006 ) detected a galaxy cluster at ffnooo = 
14 h 16 m 09 s .6, (5j2ooo = +44°44'02'.'4, coincident with the western 
shear peak, comparing archival Chandra data to optical g'r'i' 
observations in the Chandra Multiwavelength Project. They as- 
sign a redshift z = 0.427 to the cluster, designated BGV 50. In 
the same Chandra observation, Barkhouse et al.| ( |2006| l identi- 
fied another cluster, BGV 53 at a J2 ooo = 14 n 16 m 27.6 and £j 20 oo = 
+44°52'44'.'4 and a redshift of z = 0.452, which does not corre- 
spond to a bright galaxy in the Megacam image or a peak in the 
S -statistics. 

Out of the three confirmed clusters in the field, 
CL 1416+4446 not only is the only 400d X-ray cluster and 
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Fig.D.l. Like Fig. [2] but for CL 0159+0030. In the plot of rim against cnfw, a square and dashed contours denote the model 
minimising Eq.[3]if no dilution correction is assumed. 
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Fig. D.2. Like Fig. D. 1 but for CL 0230+0030. Note that, in the lower right panel, the filled circle and downward triangle, denoting 



the best parameters for the free-CNFW abd B12-models are almost coincident. 
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Fig. D.3. Like Fig. D. 1 but for CL 0809+28 1 1 . 
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Fig. D.4. Like Fig. [2] but for CL 1416+4446. Small star symbols indicate the positions of further clusters in the field, which might 
be in physical connection to CL 1416+4446 as parts of a super-cluster. 
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